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modulation  of  longitudinal  ranj'e  is  available.  .Selection  of  landinf*  site  is  possible  at  the 
expense  of  wailiny  several  orbits  until  the  earth  rotates  within  the  orbital  plane.  The  addi¬ 
tional  propulsive  weijilit  re(|uired  to  prixluce  sufficient  transverse  impulse,  which  provides 
a  lai'jte  rotation  of  the  orbital  plane,  is  impractical  svith  present  day  rocketry.  In  short, 
the  draji  mode  provides  for  a  safe  re-entry  with  little  choice  of  precise  site  selection. 

Since  t'’e  use  of  leroilynamic  drac  has  pioied  lo  Iv  a  much  more  •fficieni  means  of 
deceleraf ini'  re-entry  '.ehicles  than  propulsion,  the  invest i|iai ion  ot  the  use  of  aeroilynamic 
lift  io.icliieve  l.itcral  maneuverability  follows  naturally.  I.iftinii  re-entry  confifiurat ions 
also  offer  the  potential  of  providim:  horizontal  landinu,  which  would  furnish  better  opera¬ 
tional  capability  Ivctiu.se  lartie  recovery  task  forces  and  the  ixtrils  of  parachute  landings 
in  water  would  Ik-  elitninated.  Lifting  configurations  are  more  complex  liccause  of  the 
additional  flight -control  system,  variable  stability,  higher  rerjuired  performance,  and 
longer  flight  vlur  ition  in  the  .sensible  unvisphere.  The  duration  «>f  the  lifting  re-entry 
lie.it ing  pulse  i.s  me.isurevl  iii  lnuir.-'  Iii.-.te.id  of  minutes  as  in  the  drag  mode;  therefore, 
.iblation  technifiue.-i  are  impractical.  Instead,  the  radiatitm  ctxding  concept  i.s  used.  The 
surfai  e  is  permitted  to  reach  a  sufficiently  liigh  temperature  so  that  the  radiation  heat 
efflux  balances  the  convective  heat  influx  and  an  eiK-rgy  equi librium  occurs.  The  resulting 
high  temperautres  ret|uirc  heavy  outer  panels  to  protect  the  interior  structure.  Thus, 
the  lifting  re-cut ry  vehicle  incurs  a  greater  weight  penalty  than  the  drag  capsule  but 
achieves  the  greater  performance  of  lateral  maneuverability,  range  modulation,  and 
liori/ontal  landing.  In  short,  for  missions  requiring  lateral  maneuverability  and  precise 
site  selection  as  well  as  hori'/onial  landing,  the  lifting  mixle  is  desirable. 


Manuscript  released  !iv  author  «»n  dl  llecemlK’r  iuf(2  for  publication  as  an  .\SI)  Technical 
Documentary  Ke|Tort. 
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In  V  HI, A-  111  iIk  iMiMU'i'  Ilf  till"  lifiiiin  ili'viif,  ilic  iiiri'i'iifi  ik.’Kipnor  is  forccvl  ;«> 

I'l'i  inii.'i'  ilif  l  oiii ii',ur.ii  imi.  t  Ik-  nf  iliis  invL-st  iyiiit Inn  i.s  In  iKM'forin  ihf  nptlmiza- 

ii'Mi  ni  liliin  '  I'l -I'liirv  I'nnliv'.iiiMt inn.s  in  ,u•hk•VL■  iiia\iiiuiin  (KTfnniKinoo  fnr  ihc  minimum 

lUvinri.-  ,UK-nnii  mil  ilii-  ik-:-iymnf  iIk‘  wlik'ln.  ^vi.-  must  further  ilelincatu  tltu  design 
nlin'i  l  ive>; 

I.  \i’lik'\i.'  suiwasful  m.iojuj  iv-eivry  frnm  ciri.'iil;ir  speeds 

J.  I’rnvidi'  l.iiY'e  iaieral  maneuverahility  ami  Inngiiuilinal  range  moilulutinn 

k  I’mvivle  sat  isfaeiorv  lamlinii  charaeieristles 

•1.  I’rnteet  the  vehk'le  frnm  tile  aeriKlynamif  heating 

I'l'iivUle  llie  vehiele  svliii  ader|iiate  stahilily  In  produce  satisfactory  hamlling 
f|ualii  ies. 


Asn-  rni{-(.2-i  111') 


si  ;{  rioN  2 
constkainis 


U  iih  the  est;thlislimeni  of  the  ohjcviives  iiiul  rff|iiirfineni.s,  the  formulation  of  the 
prohlem  anil  the  conversion  of  the  physical  situation  into  the  mathematical  follow.  The 
ol)jectives  staicil  in  tlie  IntriKluct ion  will  now  Iv  separately  convert eti  into  definite 
matlictitat  ical  expressions. 


I.  Achieve  successful  manned  re-entry  from  circular  speeds 

Circular  speetl  is  2.S/»:{(i  feet  per  second.  Tor  liftinf;  re-entry,  the 
following:  e<(uilil>rium  jilide  expression  applies  (liquation  136,  Appendix  I): 


(1) 


A  suhtle  implication  in  this  statement  is  that  man  and  vehicle  have  been 
inserted  into  a  near-earth  circular  orbit  by  an  existing  or  programmed 
Itooster.  Since  rocket  Ixxtsters  are  inefficient  propulsion  tievices  (approxi¬ 
mately  MHI  poumis  of  pro|>ulsion  is  requireil  t«>  place  I  pound  into  orbit) 
and  permit  small  limiteil  payloiids,  the  re-entry  vehicle  is  weight- 
'  limited.  Hence,  tin?  specified  weight  of  the  re-entry  vehicle  Is  constrained. 
Therefore, 

w  -  filed  volute  >  (2) 


■J.  Provide  large  maneuverability  and  longitudintil  range  modulation 

As  pointed  tnit  in  Appendix  I,  maximum  lateral  ami  longitudintil  range  may  lx? 
achieved  by  optimizing  tlw  hypersonic  l./l)  (lift-to-drag  ratio)  of  the  re-entry 
spacecraft.  Thus,  tlw  hypcr.sonic  l.'l)  should  be  maximized  within  tlx?  specified 
constraints.  I'he  hypersonic  aerodynamic  characteristics  for  a  gencralizeil 
configuration  arc  cxpres.scd  in  liquations  (IifO)  thnnjgh  (216),  Appendix  III. 

3.  Provide  .stitisfactory  lamling  characteristics 

For  a  successful  landing,  the  sinking  spe-ed,  horizontal  velocity,  and  attitude 
must  be  limited  to  rcaitontibly  low  values  at  touchilown.  Appendix  V  shows  these 
minimum  acceptable  subsonic  vehicle  charac/teristics,  which  permit  a  satisfactory 
landing  to  Ik*  expressed  as  follows; 

*  ■  s  50  pounds  p*r  squOft  foot  (3) 

S  » 

4.  Protect  the  vehicle  from  acrixJynitmic  heating 

The  most  efficient  and  practical  technique  for  protecting  a  glide  vehicle  from 
aerodynamic  heating  is  the  radiation  cotding  concept  tlescrlhed  In  Appendix  11. 

The  vehicle  must  he  constructed  and  flown  so  as  not  to  exceed  the  allowable 
material  temperatures  Ity  means  of  unacceptable  attitmle,  performance,  or 
geometry.  Thus,  the  maximum  vehicle  temperature  Is  limiteu.  Four  critical 
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loiMiions  affoctinj'  geometry  will  lx;  specified  and  may  lx;  computed  from 
i-quations  (15.S),  (164),  (170),  and  (172)  of  Appendix  II.  Thus,  nose,  wing  ami  fin 
leading  edge;:,  and  lower  surface  maximum  temperature  arc  constrained.  There¬ 
fore, 


Tn  = 

fixed  voi ue 

w 

- 

fixed  VO lu e 

(^) 

Vle  ■ 

fixed  vol u e 

(f>) 

Tls  = 

fixed  voiue* 

(7) 

5.  I'rovide  the  vehicle  with  adequate  stability  to  produce  .satisfactory  handling 
qualities 


I’ilot-ratcd  simulator  studies  have  indicated  that  the  most  desirable  handling 
qualities  are  obtained  when  the  short -period  mode  of  the  rigid  body  has  a 
frequency  of  0.7  cycle  per  second  and  a  tiamping  ratio  of  0.7  (Appendix  IV). 

In  hypersonic  flight  negligible  aerodynamic  damping  occurs,  which  necessitates 
changing  the  handling  quality  criteria.  Therefore,  very  low  frequencies  are 
required  for  the  pilot  to  damp -out  quickly  oscillations  by  “out  of  phase”  control 
modulation.  Zero  frequency  may  lx;  obtained  by  imposing  neutral  aerodynamic 
sttibility  or 

=  0 


(«) 

(9) 


Cn/ 


=  0 


(10) 


as  M  approacnes  infinity. 

Appendix  V  .shows  that  a  fin  size  of  It)  percent  of  the  reference  planform  area 
is  required  to  provide  .satisfactory  subsonic  directional  stability.  Hence, 


The  problem  becomes  one  of  optimizing  the  geometry  of  the  vehicle  to  maximize 
the  hypensonic  L/D  when  the  specified  11  constraint  equations  are  used.  A 
.solution  will  be  obtained  in  the  following  .sections  by  first  determining  the 
component  geometry  that  maximizes  hypersonic  L/H  and  then  combining  these 
components  into  a  composite  configuration  and  optimizing  the  complete  geometry 
to  satisfy  the  11  constraint  equations. 
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SMc  rioN 


OP  I'IMI/.A  I  ION  Ol'  C'OMI’ONI'N'I'S 


Since’  ;in  optimi/niion  of  the  vetiicle  as  proposed  is  so  complex,  we  shall  pursue  '^hc 
prolilem  hy  analy/ini'  confiy;urai ion  components  to  determine  maximum  performance,  then 
Ity  comliininjt  tlie  pieces,  tind  by  checking  for  compatibility. 

We  shall  attempt  to  ttiaximi/e  Itypersonic  performance  characteristics  within  the 
constraints  previously  outlined  in  the  following  order: 

1.  1  .ower  surface 

2.  I’lanform 

i.  V'ertical  fin 

4.  Nose  and  leading  edge  geometry. 

With  this  itccomplislied,  the  entire  cotnposite  configuration  will  lx;  optimized  later  in 
Section  4. 


S-A.  I  .o'.ver  Surface 

The  objective  of  this  section  is  to  determine  the  shape  of  the  exposed  lower  surface  to 
produce  Uie  maximum  hypersonic  L/1).  .\  coordinate  system  is  established  in  the  con¬ 
ventional  wind  axes  system  (I'igure  1). 


I'o  assist  us  in  providing  directional  stability,  we  establish  a  plane  of  symmetry.  The 
plane  that  contains  the  x  and  z  axes  is  chosen  as  the  plane  of  symmetry  for  which  a 
inirror  image  relationship  exists  Ixttween  surfaces  on  alternate  sides  of  the  x-z  plane. 

The  equation  for  a  surface  is 


z  -  f  (  X  ,  y  )  X  0  , 

I'he  er|uation  for  tlie  outward  normal  to  the  surface  is 

7  {  I  -  f  ) 

"  =  - T  • 

|V  (z  -f  >1 

lixpansion  of  liquation  (lA)  produces  the  following  relationship: 

dz  dz 


(12) 


(13) 


ft 


(14) 


Figure  1,  Generalized  Surface  in  Wind  Axes  Coordinate  System 


riic  cloiTicntal  area  is  expressed  as  follows: 


dA 


*  (fr)‘  '■ • 


dr  .* 

y 


(15) 


As  indicated  in  Appendix  111,  only  tne  continuum  regime  requires  investigation  so  that 
modified  Newtonian  floy  may  be  used  in  ascertaining  the  pressure  on  the  surface  in 
hypersonic  flight  as  teifresentcd  by  the  following  equation: 

Cp  :  k  cos*  9  (16) 

where  9  is  the  angle  Ixitween  the  surface  normal  and  the  velocity  vector,  or 

V 

,  •  n  =  cos  9  .  (^0 

I  V  I  - 

Therefore,  by  use  of  liquation  (14)  and  since 


V 


(18) 
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and  n  =  constant  for  given  flight  conditions. 
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F  =  C. 


Hz/  (  I  -  X,  z,  ) 


I  +  j/ 


(35) 


Rc|u;uion  (35)  is  ii  cubic  eciuation  in  with  constant  coefficients.  Thus, 


*  ^4  • 


(36) 


anti  since  /  -  0, 

V 


z  :  *  +  C, 


(37) 


I  iicrefore,  the  surface  that  produces  a  relative  maximum  hypersonic  L/D  is  a  plane 
parallel  to  the  y-axis  at  an  attitude  of  a  such  that  =  -tan  o. 

This  particular  solution  remains  to  be  proved  as  the  one  desired.  The  solution  was 
obtained  by  inspection  of  Rqua  ion  (27),  i.c.,  F,^  =  C^,  Equation  (31).  It  followed  that 

X 

Zy  =  0  for  this  situation.  Therefore,  a  maximum  extremal  for  the  following  equation  was 
obtained: 


-// 


(  I  -  X,  )  dn  dy 

rrr? 


(38) 


The  more  general  equation  for  which  a  solution  is  desired  is  the  following  *. 


=// 

A 


k  z/  (  I  —  X,  z,  )  (J»  dy 
1+  z/  +  z/ 


(39) 


But  Zy  is  real  in  the  physical  problem;  therefore,!/  >  0.  Since  r/ appears  in  the 
denominator  of  Equation  (39),  it  is  apparent  that 


r,  <  Iq  for  7^  ^  0  (^0) 

and 

I,  =  Xg  for  Zy  =  0  .  (41) 


Therefore  since  the  maximum  value  of  l^  is  desired,  and  the  solution  of  produces 

the  maximum  extremal  of  I, ,  the  solution  obtained  by  inspection  (F  =  C, ,  Equation  31) 

*  \  ^ 

is  the  desired  solution. 
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I’lanivxiii 


Hy  iisiny;  ilic  results  of  the  optiini/.eii  surface  (i.e.  flat  Ixittom  surface),  we  shall 
uptinii/e  the  planform  peontetry.  Hccause  of  the  heating  constraint,  the  edges  of  the 
planforni  must  possess  a  finite  radius  to  niaint  in  specified  nose  and  leading  edge 
allowahle  temperatures  (sec  liquations  169  and  170  of  Appendix  11).  These  radii  should 
Ik;  as  small  as  tlic  heating  constraint  permits  to  produce  minimum  drag.  The  expression 
for  the  nose  and  leading  edge  aeroilynamic  coefficients  (from  Appendix  Ill)  Is  as  follows; 


and 


>«  tt 
2S 


(42) 


The  ■•I'trihurion  to  the  total  lift  coefficient  is  negligible  for  the  situation  in  which  the 
nose  and  leading  edge  areas  are  small  as  compared  to  the  total  planform  area.  Since  the 
nose  and  leading  edges  arc  temperature-limited  (see  constraint  Rquations  4  and  5), 

Rle  =  Ro  A  (43) 

wlierc  11^  =  constant  for  specified  flight  conditions  and  temperature  limit.  From  the 
planform  geometry 

2  A  cot  A  t  U  . 

coi  A  -  - )  for  lorge  A  volues  .  (44) 

14-  cot*  A  '  ' 


Thus  when  liquations  (42),  (43),  and  (44)  are  combined, 


>»  ”• 
2S 


0k  Rq 
”  3S 


/  (-M 


dK 


(45) 


where 


S  = 


(46) 


The  L/D  of  the  configuration  can  be  maximized  by  minimizing  for  prescribed 


o 

planform  surface  area,  S,  ard  fixed  nose  and  leading  edge  temperature.  The  integral, 
which  must  be  rendered  an  extremal,  is  as  follows: 


=  Cq  +  ^2  S 
0 

where  from  Equations  (45)  and  (46) 


/  G  dx 
0 


u 


G 


(47) 


(48) 


1C 
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riio  cnil  cc'iulitions  tl'-t  pertain  are  as  follows: 

at  the  nose, 


X  -  c 
y  =  0 

at  the  trailing  edge, 

»  =  0 


xe  d  end 


free  end 


liuler's  ec|uation  for  G  =  G(x,y,y')  is 


or  when  IJlquation  (48)  is  used. 


— -  6„  *  -  Gy 

d  X  y  y 


—  (5y-)  =  X,  . 


(49) 


(50) 


Integrating  with  respect  to  x  and  using  the  free-end  boundary  condition  (G  ',  =  0  at  x  =0) 
gives  y 


ay  4 

—  -  (  0. 2  X,  X  )  • 

dx  2 


(51) 


Integrating  again  with  respect  to  x  and  using  the  fixed-end  boundary  condition  (y  “  0, 
X  =  c)  produces  the  following: 


y 


{ 0,2 


^2 


-L 

4 


) 


(52) 


X^  may  now  be  determined  from  the  condition  that  S  must  be  constant  when  Equation  (46) 

is  integrated.  Thus  the  planform  for  minimum  is 

o 


y 


0.9S 

c 


(53) 


However,  y  =-^  at  x  =  0  and  S  =-^  .  Hence, 

4  1*0 


(54) 


Thus,  for  fixed  area  and  a  constant-temperature  leading  edge,  the  minimum  drag 
planform  is  a  5/4  power  curve.  The  resulting  configuration  is  not  extrem.ely  practical 
from  a  fabrication  standpoint  since  both  a  variable  leading-edge  radius  and  variable 
sweep  are  involved.  In  addition,  the  combination  of  fin  and  wing  together  at  angle  of 
attack  must  lx;  considered.  Because  of  the  Incompleteness  of  this  analysis,  the  5/4 
power  curve  should  be  used  only  as  a  guide  to  enable  the  optimization  of  the  composite 
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confii'uration.  I'or  example,  a  Bimpie  tieira  may  lx:  used  to  approximate  the  5/4  power 
curve  anti  does  not  require  variable  leading-edge  radii  or  sweep  (note  Figure  2), 


Figure  2.  Geometry  of  Minimum  Drag  Planform 


3-C.  Taper  Ratio 

When  the  stability  constraint  is  imposed,  it  becomes  apparent  the  profile  area  aft  of  the 
center  of  gravity  (c,g.)  is  required  to  provide  directional  stability  at  all  speeds.  Conven¬ 
tionally,  vertical  fins  in  the  most  aft  position  have  been  used  to  provide  this  stability.  A  ven¬ 
tral  fin  would  most  efficiently  provide  the  desired  stability;  however,  the  heating  is  so  severe 
that  the  dorsal  position  is  the  only  presently  acceptable  location.  Dorsals  also  provide 
lower  drag  than  a  ventral  (for  the  same  material-temperature  limit)  because  of  the  greater 
effective  sweep  when  operating  at  angle  of  attack  greater  than  zero.  Since  the  upper  sur¬ 
face  is  shielded  hypersonically.  two  dorsal  fins  located  outboard  in  the  most  aft  position 
are  required.  As  previously  shown,  the  planform  of  fixed  area,  temperature,  and  chord 
that  presents  the  minimum  leading  edge  drag  is  close  to  a  delta  shape  (Figure  2).  For  this 
reason  the  vertical  fin  will  be  configured  in  this  optimum  delta  shape.  Thus  the  optimum 
configuration  possesses  a  flat  bottom  (Section  3-A),  a  delta  shaped  planform  (Section  3-B), 
and  two  delta  dorsal  vertical  outboard  fins.  The  vertical  fins  can  be  attached  effectively, 
however,  by  clipping  the  planform  delta  as  shown  in  Figure  3.  Note  the  root  chord  of 
length  c,  taper  ratio  m,  wing  sweep  A>  and  fin  leading  edge  sweep  Ap.  This  geometry  will 

be  optimized  for  maximum  L/D,  by  imposing  the  constraints  established  in  Section  2. 
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FlRurc  3,  Geometry  Model  for  Optimizing  Taper  Ratio 
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Tho  axial  force  coefficient  of  all  componeiit.s  for  the  descrilied  configuration  is  as 
follows  (see  Appendix  111): 


TT  k  Rk 


M  jkRirb  2 

-  (  I  +  cos  a  )  + - = —  tcosA.  +cosA  cos  a) 


'A  - 


8S 


3S 


Sr  1  2  8k  R IT  he  2  * 

+  2k  — cos  a  +  - -  cos  (Ap  +  otl 


,,  0.45  cos  a  4.65  1 

8  Jw 

f  '' 

z,z 

1  sin  a  cos  o 

1 

.  '  » 

f 

1  10,000  t 

s 

>  r  .0» 

(55) 


Frttn.  the  heating  ec|uarions  of  Appendix  11 

W  ^  8.  025 


‘N 


SC, 


'N 


(56) 


w 


sc. 


,  0.675  X  t. 

(1 - ) 

\  «LC  ' 


(57) 


where 

cos*  Ag  •  l-sin*Acos*a  for  wing  leading  edge 
cos^Ag  =  cos*  (  Ap  +  a)  for  fin  leading  edge  . 

The  critical  design  point  for  heating  is 

=  21,170  feet  per  second  (Appendix  11). 

When  the  equilibrium  glide  equation  (Appendix  1)  is  used,  the  Reynolds  number  may  be 
ascertained. 


V  c 


2W  c 


s  Cl 

At  critical  heating  and  for  «3  x  lO'"  lb-scc./ft2 

V  c  105  w  c 
'  s  Cl 


(58) 


(59) 


From  the  geometry  of  Figure  3,  the  following  quantities  may  be  obtained  (for  small  X 
values): 


(  I  4-  rn  )  , 


(60) 
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where 


S 

jlf.  .  /  14  m  \  Sp 
b  '  m  '  s 


(61) 


arnl 


cot  A  = 

4  V  I  —  m  ) 


(62) 


cot 


^  lE.  (_L±I!L) 

Z  S  '  m  '  • 


(63) 


Siiii^^iitution  of  Ixiuations  (36),  (37),  (34),  (60),  and  (61)  into  liquation  (55)  results  in  the 
following: 

TT  k  ( I  +-  cos  a  I*  /  vw  \*  /  8.025  x* 


8tV  C| 


(1)  {^) 


kb  ,  Wf  ,5  675  r  ,  ,  .  ,  z 

~  \~5 - /  i®*  Ag  (cos  A,  4  cos  A  cos  a  ) 


3W  C, 


SSp  /  l  +  m\  4,  1  Sc,  - 

—  [— — )  cos  (  Ap  4-  a  )  j  +  2k  X*  cos*  a 


(64) 


4-  Y  “Y  (“iiovv')  °  “  cos  *a  I  - - -  • 

When  liquation  (64)  is  written  in  functional  form  with  the  aid  of  liquations  (62)  and  (63), 


Ca 

NOSE 

(  w  . 

w 

s  ’ 

Q  k  Gt 

) 

ca 

LE 

(  b  , 

w 

T 

-  .  Ole  ' 

Sf  , 

,  a,  m  ) 

(65) 

Ca 

FIN 

s 

.  X. 

a  ) 

c, 

/  Sw 
'  s 

.  b. 

W,  a  ,  m  ) 
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I lowovt'r,  fniiii  tlio  constraints  iniposcJ  its  nivcn  in  Section  2 

constroini  equation 


w 

= 

con  s  f  a  n  f 

(2) 

fc 

—■  ..1  ...1  r  constont 

(31 

>/50 

const  am 

(4) 

= 

consto  n  t 

(5,6  ) 

2.26  =  constont 

(7) 

sf 

0.  1  • 

(II) 

=  2,  o  =  and  X  determined  from  constraint  lifiuiition  (*?)  then  liquation  (65) 

=  ('  ^(m).  To  maximize  L/O,  one  should  optimize  the  taper  ratio  (m)  to 

minimize  Only  the  second  and  fourth  terms  (leading  edge  and  skin  friction  terms)  of 

liquations  (6-1)  and  (65)  depend  upon  m.  However,  the  skin  friction  expression  depends 
only  mildly  upon  m  as  follows; 


,,  W 
I- or 

liocomes  C 


r . . . — . 

m 

1  -  m* 

0 

1  .  000 

0 . 2 

1  ,  039 

0 . 4 

1  .  052 

0  6 

1 . 058 

0 .  e 

1 . 060 

1 . 0 

1.061 

Therefore  for  constants  W,  b,  and  a,  the  laminar  skin  friction  coefficient  is  virtually 
independent  of  taper  ratio  (only  ±3  percent  variation).  Hence  the  optimum  taper  ratio 
must  be  determined  from  the  condition  for  minimum  leading  edge  drag. 

The  leading  edge  drag  component  of  Equation  (64)  may  be  written  as 


( —55 — ^  Co,  c-  ■  ico»  A,  +  cos  A  cos  a  1* 

'kb  Ra  ' 


esr 


(-L^)  co.^Ap4.a, 
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where 


w 


5.  675 


(67) 


R 


o 


SCl 


( 


‘•le 


When  the  geometric  relationships  (Equations  62  and  63)  are  used, Equation  (66)  may  be 

solved  for  various  values  of  m.  =  0.1,  and  for  fixed  values  of  and  a.  In  the  region 

of  interest  (0.5  2  and  10“  <  n  <  25“),  numerical  calculations  show  that  a  taper  ratio 

between  0.30  and  0.35  produces  minimum  total  leading  edge  drag  (Figure  4).  Therefore, 
an  average  value  of  0.32  for  optimum  taper  ratio  was  selected.  Since  the  c.g.  must  be 
located  near  the  area  centroid  of  the  lower  surface  (about  0.35  c  from  the  trailing  edge) 
to  satisfy  the  pitch  stability  constraint  (Equation  9),  the  optimum  taper  ratio  of  0.32 
advantageously  places  the  entire  vertical  fin  aft  of  the  c.g. 

">ji-  the  configurations  to  be  investigated  herein,  a  taper  ratio  of  0,32  was  adopted. 

Sp 

However,  note  that  for  values  of  ,  and  a,  which  are  much  different  than  those 

stated  above,  sizable  differences  in  optimum  taper  ratio  values  will  result. 


3-D.  Nose  and  I.eading  Edge 

The  flight  corridor  of  re-entry  vehicles  is  limited  by  its  ability  to  absorb  or  reject 
the  re-entry  heat  load.  The  heat  pulse  applied  to  a  ballistic  vehicle  during  re-entry  is 
of  such  short  duration  that  cooling  by  ablation  allows  the  vehicle  to  survive.  Because  of 
long  re-entry  times,  however,  the  only  presently  practical  way  to  cool  lifting  vehicles 
is  by  radiating  their  heat  loads  to  space.  The  long  re-entry  time  and  small  heat  capacity 
of  a  lifting  vehicle’s  skin  results  in  equilibrium  between  heat  convected  to  and  heat 
radiated  away  from  the  surface. 

The  objective  of  this  section  is  to  determine  the  detailed  shape  of  the  nose  and  leading 
edge  that  will  minimize  the  peak  temperature.  When  the  peak  temperature  is  minimized, 
smaller  radii  configurations  may  be  employed,  which  result  in  lower  total  drag. 

As  shown  in  Appendix  II,  the  highest  heating  occurs  at  the  stagnation  region  in  laminar 
flow  at  a  velocity  of  21,170  feet  per  second  during  glide  re-entry.  Consider  the  total 
heating  rate  into  a  surface  in  radiation  equilibrium  flying  at  this  peak  heating  point 

6  =  /  q  dA  z  f  a-  J  J*  dA 

where  constant  emissivity  has  been  assumed.  An  average  temperature  may  be  defined 
as  follows: 

or 


(70) 
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Figure  4.  Wing  and  Fin  Leading-Edge  Drag  Ver«u»»  Taper  Ratio 


If  a  region  of  A  exists  for  which  T  <  Equation  (3)  shows  tliat  another  region  must 

exist  for  which  T  >  T^yg-  Thus,  the  surface  characterized  by  A,  Q,  and  c ,  which  has  a 

constant  surface  temperature  T  ■  T^^g,  has  the  minimum  peak  temperature  of  all 

surfaces  characterized  by  A,  Q,  and  c.  Since  q  Is  related  monotonlcally  to  T.  a  consent 
T  ■  T^yg  implies  a  constant  q  ■  Therefore,  for  a  surface  characterized  by  A,  Q, 

and  r,  which  Is  maintaining  thermal  equilibrium  by  radiation,  the  minimum  peak  tempera- 

4 

ture  will  occur  when  the  local  heating  rate  on  A  is  a  constant  q  ■  co  T^yg*  ITius,  the  prob¬ 
lem  of  finding  the  optimum  nose  and  leading  edge  geometry  is  reduced  to  that  of  finding 
a  surface  of  constant  heating  rate. 
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I.ees  (Heferenco  2)  dcBcribcs  the  laminar  convective  heat-transfer  rate  near  the 
stagnation  point  of  a  blunt  body  in  hypersonic  flight  by 


where 


F  (s) 


V,  ds  4  ^ 


F(») 


(71) 


(72) 


where  n  =  0  for  a  planar  body  (leading  edge)  and  n  =  1  for  a  body  of  revolution  (nose). 

The  static  pressure  ratio  will  be  represented  by  the  Newtonian  flow  approximation 

=  cos*  e  (73) 


with  0  as  the  angle  between  the  surface  normal  and  the  velocity  vector.  Measurements  by 
M.  O.  Creager  (Reference  3)  and  others  show  that  the  total  pressure  is  nearly  constant 
along  invisciu  oucaialinco  tciicv.ir.j,  a  biem  IvJy  •  ontour.  When  this  is  assumed,  the 
local  velocity  at  the  outer  edge  of  the  boundary  layer  may  be  related  to  the  local  pressure 
by  the  following  isentropic  relationship: 


For  hypersonic  velocities: 


(74) 


ta  I  • 


Substituting  Equations  (73)  and  (75)  in  Equation  (74)  gives 


vg 


V. 


e 


where 


(75) 


(76) 


= 


y-i 

y 
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The  velocity  gradient  may  now  be  determined 
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dv 

da 


/^/3*  [l-< 


1  2 


S  ion  d 


dg 

ds 


(77) 


The  value  ot  the  velocity  gradient  at  the  stagnation  point  is  required  for  Lees’  heating 
rate  equation. 


where 


/  R  dVg  V 

s 

tan  6  ] 

'  dv  J 

1  -  cos  g  ' 

^  J 

ds 

d? 


=  R. 


(78) 


(79) 


When  Equations  (72),  (73),  (76),  and  (78)  are  used  in  Equation  (71),  the  hearing  rate 
distribution  may  be  obtained  for  an  arbitrary  body.  Following  Lees  (Reference  2),  the 

^  6 

assumption  is  made  that  —  =  1  in  hypersonic  flight. 


Wagner,  Pine,  and  Henderson  (Reference  4)  have  shown  experimentally  that  Lees’ 
equation  holds  well  for  a  series  of  blunt-nosed  bodies  with  noncircular  cross-sections. 

Since  Rg  appears  only  as  a  scaling  factor  in  Equation  (80),  it  may  be  assumed  unity 

without  loss  of  generality.  The  heating  rate  ratio  of  Equation  (80)  must  be  unity  to 
produce  a  constant  heating  rate  over  the  surface. 

Defining 

jl_ 

g(g)  =cos*0  (i  — co»^^  9  )  (®^) 

gives  the  equation  for  a  constant  heating  rate  surface  (after  squaring) 

g(g)  d*  »  g*(g)  y*"  (82) 

o 

Differentlaticri  produces 

JL  =  fl  2"  -  -^2-  (83) 

dg  dv  ^  y  ds 
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Thus 


_d»_ 

dd 


dg 

de 


0  2*' — ^  cos  8 


The  geometric  relationships  of  x.  y.  s.  0  are  as  follows: 

d»  ds 

— ^  =  sin  S  — — 
d&  dd 


dy  ds 

— —  s  cos  0  - — 

dd  dd 

or  ;  .n  r  substitution  of  liquation  (81)  into  liquation  (84) 


sin  8  cos  8 


d  s 

dfl' 


«  2  4-/3*  )  cos^^  e-2 
[  (.-co.2^*0)? 


/9  2**  —  —  cos*0  {  l-cos  8  )  ' 


(84) 


(85) 

(86) 


(87) 


By  simultaneous  integration  of  Equations  (85),  (86),  and  (87)  with  the  initial  boundary 
conditions  that  x  =  0  and  y  ■  0  at  0  »  0,  numerical  solutions  were  obtained  for  two  and  three 
dimensional  cases  (n  =  0  and  1,  respectively)  and  for  y  of  1.4,  Numerical  results  are 
tabulated  in  Table  1. 


In  Table  1,  note  the  small  percentage  difference  between  the  optimum  geometry  Tor 
two  and  three  dimensional  surfaces.  The  heating  rate  distribution  on  the  optimum  surfaces 
are  within  3  percent  of  the  distribution  for  a  2:1  ellipse  (minor  axis  parallel  to  the  velocity 
vector).  Therefore,  for  simplification,  the  nose  and  leading  edge  optimum  geometry  will 
be  constructed  as  a  2:1  ellipse.  From  Table  1,  also  note  that  beyond  a  slope  of  36"  the 
radius  of  curvature  becomes  negative  (a  cusp  occurs)  and  the  optimum  surface  can  no 
longer  be  maintained.  It  can  be  shown  from  Equation  (80)  that  attaching  the  optimum 
surface  to  the  remainder  of  the  vehicle  with  a  smooth  curve  of  positive  radius  of 
curvature  will  insure  that  at  no  point  will  the  stagnation  region  temperature  be  exceeded. 
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TABLE  I 

C0NSTAN1  TEMPERATURE  SURFACE  GEOMETRV 


y«l.4,  )9»  0  ,53452 


"  -  0 

n  «  1 

o 

( DEGREES ) 

X 

R* 

y 

R» 

±  iL 

R*  d®. 

X 

Ri 

y 

Rf 

_L  ^ 

R*  d0 

0 

0 

0 

1 00000 

0 

0 

(.00000 

2 

60794** 

34779*' 

99659 

.60795** 

.34854*' 

99608 

4 

.24182** 

.69342  * ' 

.98570 

.24(69** 

.69393*' 

.98438 

6 

.53893** 

10332 

9703  ( 

.53820** 

.(033( 

.96502 

e 

94527'* 

.  1 3644 

.9463  ( 

.94298** 

.  ( 3630 

.938(7 

10 

.14514*' 

.16842 

.91552 

.(4459*' 

.(6807 

.904(3 

12 

.20456*' 

.19900 

.87820 

.20345*' 

.(9837 

.86325 

14 

.27137*' 

.22796 

.83464 

.26936*' 

.22695 

.8(593 

16 

.34397  *' 

25508 

78520 

.34066*' 

.25358 

.76269 

IS 

42058*' 

.28015 

.73028 

.41549*' 

.27807 

.70409 

20 

.49930*' 

.30303 

.67032 

.49(88*' 

.30028 

.64075 

22 

.57808*' 

.32358 

.60583 

.56779*' 

.32007 

.57339 

24 

.65484*' 

.34168 

53733 

,64(14*' 

.33737 

50274 

26 

.72744*' 

.35726 

46538 

.70987*' 

.352(3 

.42963 

28 

.79374  *' 

.37029 

.39059 

.77(98'' 

.36433 

.35489 

30 

.85165  *' 

.38076 

.3(358 

.82558*' 

.37402 

.27942 

32 

.899(4*' 

.38868 

23500 

.86897*' 

.38(26 

.204(3 

34 

.93431  " 

.394(1 

(5552 

.90064*' 

.386(5 

.(2995 

36 

.95541  *' 

.397(3 

.758(4*' 

.9(933*' 

.38883 

.57809*' 

37.719 

38 

96085  *' 

.39787 

-.(2637*' 

.92429*' 

.38950 

-.(8802‘* 
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Si'CTlON  4 

Ol’  I'lMI/ATION  or  COMl’OSITI-:  CONFIGURATION 


4-A.  Introduction 

In  the  preceding  Section  3,  we  showed  that  the  maximum  hypersonic  L/D  may  be 
achieved  hy  a  vehicle  with  the  following  characteristics: 

1.  Flat  bottom  surface 

2.  Clipped  delta  planform  with  a  0.32  taper  ratio 

3.  Outboard  dorsal  delta  fins 

4.  Flliptical  (2:1)  nose  and  leading  edges 

In  this  section  the  composite  configuration  will  be  optimized  for  maximum  hypersonic 
U/0.  which  satisfies  the  11  constraint  conditions  delineated  in  Section  2. 

These  constraint  equations  are  summarized  as  follows: 


P  V* 

JL  JiL) 

SCl  '  Vg*  > 

(88) 

w 

fixed  value 

(89) 

w 

SfJ' 

r 

50  pounds  per  square  foot 

(90) 

Tn 

= 

fixed  value 

(91) 

r 

fixed  volue 

(92) 

Vle 

z 

fixed  volue 

(93) 

\s 

= 

fixed  value 

(94) 

= 

0 

(95) 

r 

0 

(96) 

C„^ 

z 

0 

(97) 

Sf 

s 

z 

0.  1  . 

(98) 

The  available  equations  from  which  a  solution  may  be  attempted  will  now  be  reviewed. 
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•J-B.  A v;i il;il>U'  llriiintions 

In  Appcnilix  11.  the  temperature  liquations  (155),  (IA‘)),  (170),  and  (172)  were  derived. 
They  are  summari/ed  here  in  functional  form: 


'  ‘-r 

.  Cl 

.  ) 

(99) 

^LE  = 

f  (f 

«  Cl 

,  ,  a  ,  a  ) 

(100) 

^FLE  - 

f{f 

.  Cl 

,  Rp  ,  Ap  ,  a  ) 

(101) 

1 

1 . 

(102) 

The  hypersonic  aerodynamic  equations  are  compiled  in  Appendix  III,  liquations  (190) 
rhrr.i;i;h  (216). 

[  •  ‘'N  •  ]  ' 


f  [Rn  • 

Rle  • 

Rp  .  , 

,  «lE  •  *l  •  "f  •  "fle  ■ 

>'LE  •  ''F  •  *  N  •  *LE  •  *L 

-^LE 

.Sp  ,Sl 

,  X  ,  A  ,  Ap  ,  s , 

b  ,  c  ,  c  ,  a  ,  V,  . 

The  aerodynamic  characteristics  as  shown  are  a  function  of  30  variables,  and  since 

only  1 1  cnnetrflinr  eqimfipns  nrp  aviillnhlp,  ohvlniisly  iuHlirlonal  mixllinry  relnrlonshlps 
iii'c  I'l'fiiili'tAl,  Pill'  lilt'  jtoomeify  nIiowh  In  I'ljiuni'  .1,  llio  followlnn  ncoineiHo  fclatlonNliliid 
may  lx;  obtained: 


t 

c 

=  s[<l-m*)  cot  A  +  m*  ton  xj 

(104) 

b 

-  2c  [ll  — m  )  col  A  +  >n  ton  X  1 

(105) 

cot  Ap 

2Sp 
t  «  * 

(106) 

m  c 

c 

2 

(107) 

3  ' 

Sw 

=  2S 

(108) 

Sl 

=  S  -  2-^lE  '’lE 

(109) 

■^LE 

b 

2  cos  A 

(110) 

Jr 

,  _ »!!£ _ 

(111) 

sin  Ap 
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''le 

r 

b 

4 

“  '’le 

(112) 

I 

b 

2 

(113) 

: 

*cg  ■ 

-  Rn 

(114) 

*LE 

= 

*cg 

-  Rle 

(115) 

= 

*cg 

(116) 

*cq 

- ^  me  cot  Ap 

(117) 

^FLE 

= 

*cg 

-  me  cot  Ap 

(118) 

»N 

= 

*cg 

(119) 

*LE 

= 

*cg 

-  ton  A 

4 

(120) 

*  F 

*cg 

-d  -  c 

(121) 

*FLE 

* 

*cq 

-  ( 1  —  m  )  c 

(122) 

_ 

c 


*cq 


-  1  + 


m*  ton  X 
2(1 -m)  cot  A 


I  + 


m*  ton  X 
( I  -  m  )  cot  A 


(123) 


In  addition,  the  kinematic  viscosity,  v,,,  may  be  expressed  as  a  function  of  altitude 
as  follows; 


®  5  xlO*^  lb  soc./tt*  in  tn»  altitude  roigon  of  intorast 


or 

=  t  (  /»^  ) .  (124) 

Therefore,  when  the  geometric  relationships  are  used,  the  11  constraint  equations  and 
the  equation  for  may  be  expressed  by  the  following  independent  variables: 


1. 

W 

5. 

Rf 

9 

m 

13. 

2. 

W 

S 

6. 

A 

10. 

*cg 

14. 

a 

3. 

«N 

7. 

Ap 

II. 

zcg 

4. 

Rle 

8. 

\ 

12. 

V 

V 
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Since'  12  criiKUions  exist  with  14  varlaliles,  two  aJUitional  ef|uatlon8  or  conditions  are 
refill i red.  The  velocity  for  which  tlie  peak  aerodynamic  heating  occurs  (V  =  21,170  feet  per 
second,  litiuation  16!5  of  Appendix  II)  is  the  critical  design  point  and,  hence,  the  vehicle 
must  he  constructed  so  as  not  to  exceed  the  temperature  at  that  point.  Thus 

V  =  21,  170  feat  per  eecond  (1^^) 


is  one  additional  condition. 


From  the  preceding  section,  the  optimum  taper  ratio  (m)  for  minimum  leading  edge 
drag  was  ascertained  as  0.,'12.  lienee,  the  final  condition  is  that 


m  =  0.  32  . 


(126) 


1  he  problem  now  is  properly  set  whereby  the  11  constraint  equations  and  the  maximum 
L/i)  cc'iidition  are  expressible  by  12  independent  variables.  When  these  11  equations 
arr'  •lived  simultaneously  and  the  optimization  process  is  imposed,  a  solution  may  be 
obtained. 


4-(!'.  Numerical  Solutions 

llecause  of  the  complex  interrelationship  of  the  11  constraint  equations,  an  iteration 
procedure  was  established  to  obtain  a  solution  of  the  system.  The  aerodynamic  equations 
were  programmed  on  the  IBM  7090  computer  at  ASD,  The  geometric  quantities  were 
computed  for  the  condition  of  fixed  weight  and  temperature  as  shown  in  Figure  5  in 
block  diagram  form  (trial  values  of  x^^,  X,  and  were  required  to  initiate  the 

program).  A  complete  set  of  aerodynamic  characteristics  was  obtained  from  the  IBM 
7090  computer  for  every  degree  of  o  from  7*  to  30*.  By  inspection  of  the  readout  dat  i, 
the  angle  of  attack  for  was  ascertained.  The  new  radii  were  computed  for 

the  fixed  temperatures  and  for  the  adjusted  C,  and  “q-j  values  (constraint  Equations 

opt  ^ 

‘/I,  92,  93).  A  new  equilibrium  altitude  was  computed  from  constraint  Equation  (88).  The 
vehicle  was  rendered  neutrally  stable  (con.srr.iint  Equations  96,  97,  95)  by  relocating  the 
center  of  gravity  and  adjusting  the  fin  toe-in  angle  (X)  as  follows: 


^  *cq  _ 

'  ^n7 

Cn« 

AX  r - r - 

Sp 

8  — ^  — —  cos  a 
b  S 

and 


(127) 


(128) 


(129) 
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Figure  5,  Block  Diagram  of  Configurational  Optimization  Process 


Additional  geometric  quantities  were  adjusted  to  conform  with  the  four  constraints  of 


0.1,  Y  '  through  Tj 


50  pounds  per  square  foot  and  W  =  fixed 


(Equations  98,  94,  90,  89),  This  procedure  was  continued  until  satisfactory  convergence 
was  obtained.  The  criteria  for  satisfactory  convergence  was  that  the  aerodynamic  center 
in  all  three  axes  he  within  I  percent  of  reference  length  (c  or  b)  from  the  c.g.  and  that 
the  design  angle  of  attack  (®Qpj)  he  within  one  degree  of  the  iterated  value.  Approximately 

10  successful  iterations  were  required  to  produce  satisfactory  convergence.  It  was  also 
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ascortalncil  iluit  most  siiccL'ssful  convern'jncf  w.i.s  ohi.iincil  hy  ailjusiinn  ‘'nc  riuantity  at 
a  tiiitc  in  iho  ftillowin^  sefiuomi': 


2 •  A  *cg 
3  AX 

Solutions  were  ol)t,.!.'i‘i|  at  three  temperatures  (or  svlnp  loadinps)  for  two  vehicle  weights 
of  tO.OOO  and  100,000  pountls.  Since  four  temperature  locations  were  examined  (nose,  wing 
leading  edge,  fin  leading  edge,  and  lower  surface),  some  assumption  concerning  the 
material  construction  was  required.  To  minimize  the  number  of  cases  to  lx?  examined,  we 
assumei.i  that  the  lower  surface,  wing  leading  edge,  and  fin  leading  edge  were  constructed 
of  the  satiie  materials.  That  is 


=  n 


Lf 


T  LE 


(130) 


In  addition  we  assumed  that  the  minimum  lower  surface  temperature  as  computed  from 
the  wing  loatling  value  (liquation  172,  Appendix  11,  or  Figure  IS)  was  the  design  limit.  In 
other  words  no  excess  temperature  capability  was  built  into  the  vehicle  and  hence  it  is 
considered  to  lx?  a  minimum  temperature  design. 


The  nose  of  the  vehicle  rcpre.sents  a  small  portion  of  the  vehicle’s  weight  yet  influences 
the  value  of  (l./D)  greater  than  any  other  component  of  similar  area.  Therefore,  a 

higher  temperature  capability  for  the  nose  than  the  remainder  of  the  configuration  results 
in  huge  dividends  in  hypersonic  L/D,  with  only  a  minor  penalty  in  weight.  As  stated  in 
Appendix  II,  nose-cap  materials  have  temperature  capabilities  in  the  vicinity  of  4000*F 
while  lower  surface  materials  of  lighter  weight  arc  limited  to  about  30()0“F.  For  this 
case. 


^  ^LS  ' 

Hence,  in  all  numerical  cases  investigated,  the  nose  was  assumed  to  possess  a  heating 
rate  limit  that  was  three  times  as  great  as  the  remainder  of  the  vehicle. 


The  characteristics  of  the  optimized  vehicles  for  the  conditions  cited  are  shown  in 

Tables  2  and  3.  The  values  of  (L/D)  versus  W.  for  these  two  classes  of  weight  are 

'  'max  S  “ 

presented  in  Figure  6.  Note  that  the  L/D  increases  with  weight  and  decreases  with  wing 

W 

loading  although  the  temperature  limits  were  constrained  to  increase  with  (as  dictated 


by  Equation  172,  Appendix  II),  Hence  the  decrease  of  L/D  with  increasing  wing  loading 
would  be  even  more  drastic  for  a  constant  temperature  constraint.. 


After  the  optimization  process  is  complete,  there  are  always  interesting  "tradeoffs” 
that  should  be  investigated  to  determine  the  penalty  incurred  by  various  constraints. 

For  this  reason  the  maximum  L/D  was  recomputed  for  the  optimized  configurations  with 
the  exception  that  the  following  quantities  were  removed; 


Fin  —  0;  no  stability  constraint 
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TABLE  2 

CHARACTERISTICS  OF  THE  OPTIMUM  LIFTING  CONFIGURATIONS  :  W  »  10,000  POUNDS 
m  s  0.  32,  Sp  /S  =0. 1  ,  V  s  21,170  ft  /  f«c.,  ond  b  «  14.  14  ft 


W/S  (  lb/  «q  f  t  )  at 

00 

CM 

CM 

52.6 

76 . 7 

A  (d«gr««t  ) 

82.8  1 

69.2 

57.9 

Ap  (d«4r««9) 

76.4 

55.2 

42.0 

X  (d«grt««) 

5.  90 

8.  50 

7.80 

Rn  (ft  ) 

1 . 54 

0.  67 

0.34 

Rle  (ft: 

0.51 

0.60 

0.6  1 

Rp  (ft  ) 

0.03 

0.  28 

0.29 

c  (ft) 

n 

23.1 

*  eg 
c 

0.  62  1 

0.615 

0.625 

^cg 

b 

0.  1  73 

0.  156 

0,  1  77 

R,  X  10'* 

1 .02 

0.654 

0.399 

T^e  CR) 

2750 

32  50 

3510 

T^  CR) 

3620 

4280 

4610 

Oopt  (d«gr«<s) 

14 

17 

23 

(L  /  D  )niax 

2.08 

1 . 75 

t.29 

0.  <  5P 

0.2  19 

0.369 

1 

-  0.  205 

-  0.340 

-  0.365 

Cm  (parcant  ) 

-  0.  90 

0.  27 

0.82 

Cma  (pa'^canO 

-1.56 

-  1.93 

0.00 

Cn^  7Cy^  (parctnt) 

0.65 

-0.  3  5 

-  0.68 

/Cyp  (parcant) 

0.05 

-  0.  54 

0.2  7 

20 


TABLE  3 

CHARACTERISTICS  OF  THE  OPTIMUM  LIFTING  CONFI GUR AT  IONS  I  Ws  100,000  POUNDS 
m  =  0.  32  ,  Sp  /S  =0.1  ,  V  =  21,170  fl/  sec.,  ond  b  »  44. 72  M 


W  /  S  ( 

lb  /  sq  ft  )  at 

17.6 

45.7 

7  1.0 

A  ( degrees  ) 

82.8 

69 . 2 

5  7.  8 

Ap’  (degrees) 

7  7.  0 

55,  8 

4  1.9 

^  (degrees) 

2  75 

4.  30 

4.00 

R^  (ft  ) 

3.  30 

1.14 

0.48 

Rle  ('•  > 

0.  63 

0,87 

0.  77 

Rp  (»t  ) 

0.  05 

0.  67 

0.  51 

c  (  M  ) 

21  9 

79.  Z 

49.6 

*eg 

0.6  42 

0.  622 

0.645 

"  “c 

*  CO 

0.  1  86 

0.  162 

0.224 

b 

R,  X  10' * 

•  • 

5.33 

3.03 

1  .6  2 

2610 

3160 

3450 

e'/"  Tp,  CR  ) 

3440 

4160 

4540 

Oopi  (degrees) 

10 

13 

18 

3.  26 

2. 54 

1.88 

Cn 

0.087 

0.  133 

0.235 

-  0.  075 

-  0.  1  48 

-  0.  1  55 

/Cff  (percent) 

0.  63 

-  0.  03 

0.  80 

Cma  '‘^Na 

0.  42 

-  0.  2  3 

1.2  1 

/Cy^  (percent) 

0.  32 

-  0.  2  3 

-  0.69 

/Cy^  (percent) 

0.  64 

0.2  2 

-  0.99 

30 


le 


cl 

ISI 

cl 


t  IS 

s. 


q  ft) 


Versus  Wing  Loading 


ASP-  l()L> 
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l-l).  Siirl'.ii-i- 

In  iho  pniMgraphs  iho  opiinuitn  vo'iick-  cliaractc-riKtics  for  a  lifting  re-entry 

vehk'le  were  ui'ierinineil.  TlironulMiiit  the  analysis  only  the  planform  and  lower  surface 
i-eomelry  wv-i'e  discussed.  The  lower  sui  lace  geometry  is  deslRned  by  the  hypersonic 
characterisi ics  leaving;  the  upper  surface  free  for  the  designer  to  acliievo  satisfactory 
low -speed  eharacierist ics  while  meeting  the  volume  and  Ixiostcr  attachment  constraints. 

A  c.tmhered  upper  surface  that  Utcates  the  Itwv-.speeil  center  of  pressure  aft  of  the  liyiier- 
sonic. tlly  constrained  c.g.  and  minimi/es  profile  drag,  achieves  optlmutn  low-speed 
perfortnance  and  stability.  llow<-ver,  the  volutne  and  attachment  constraints  may  alter 
this  conutur  in  various  tlegrees.  A  cttmplete  analysis  of  this  problem  Is  not  wtirranted 
here,  in  th.it  infinite  possibilities  exist;  however,  a  typical  complete  configuration  is  shown 
as  an  example. 


One  cf  the  l()0,(K)()-pound  configurations  =  4.S.7  pounds  per  square  foot)  is  shown 

in  !■  'i;r<‘  o  with  a  cambered  r.■odified  C'lark-Y  profile  upper  surface.  All  constrained 
geotnetric  c|uantitie.s  are  shown,  and  a  cylindrical  aft  fuselage  is  depicted  unexposed  to 
the  aerixlynamic  heating  to  contain  the  payload  volume  and  c.g.  while  providing  an 
interface  for  Ixtoster  attachment. 


M 


Figure  9,  Typical  Optimum  Geometry 
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l-l'..  Hanyic  MtKlnlaiion 

I  hiouj^hutu  iho  oiuire  invest ij'at ion,  ilie  veliicle  lias  lieen  Ucsij'neil  for  one  aitituiie;  that 
is,  ilie  angle  of  aiiaek  ’^”•‘*1  llL'^'ilng,  sinhility.  Internal  heat  protection,  control 

surfvioe,  sizing,  etc.  liave  lieen  determined  for  tliis  condition.  However,  L/l)  must  he 
adjusted  from  the  maxitmim  value  to  some  lower  value  to  modulate  longitudinal  range. 

1  he  degree  of  range  moilulation  depends  upon  the  ratio  of  the  maximum  l./H  to  the 
minimum  value,  l./l)  modulation  may  Iv  achieved  hy  opening  a  speed  brake  to  increase 
drag  at  ctinstant  angle  of  attack.  Symmetrical  deflection  of  the  rudders  may  accomplish 
this  feat  if  sufficient  temperature  canability  exists.  Table  4  expres.ses  the  values  of 
l./l)  for  rudder  deflections  of  4.S'’  for  the  configurations  studied.  The  rudder  area  was 
assumed  to  lx*  one-lialf  of  the  fin  area. 


table  4 

VALUES  OF  L/0  FOR  RUDDER  DEFLECTIONS 
Sp  •  0.5  Sp  and  Sp*  0.1  S 


w 

— 

W/S 

®opt 

Sr.O* 

8p.45- 

( POUNDS) 

(LB  PER  SO  FT) 

(DEGREES) 

( L/0)mox 

( L/D )min 

10,000 

22  8 

1  4 

2  080 

0.703 

1  0,000 

52.6 

17 

1.751 

0.802 

IQ, 000 

76.7 

23 

1.287 

0.805 

100,000 

17.6 

10 

3.263 

0.533 

100,000 

45.7 

13 

2.543 

0.701 

100,000 

71 .0 

18 

1.883 

0.868 

.\Sl)-I'nR-(i2-lll)2 


SI-CI'ION  5 
('ONCLDSIONS 


The  objective;  of  tliis  investi^jation  was  to  optinii/e  lifting;  re-entry  configurations^ for 
the  following  set  of  requirements: 

1.  Achieve  successful  manned  re-entry  from  circular  speeds 

2.  Provide  large  lateral  maneuverability  and  longitudinal  range  modulation 

3.  Provide  satisfactory  landing  characteristics 

t.  Protect  the  vehicle  from  the  aerodynamic  heating 

5.  Provide  the  vehicle  with  adequate  stability  to  produce  satisfactory  handling  qualities 

These  requirements  were  converted  into  the  11  property  constraints  for  the  vehicle.  With 
the  constraint  equations  formulated,  the  configuration  was  designed  for  maximum  hyper¬ 
sonic  L/D  after  it  was  shown  that  the  greatest  lateral  and  longitudinal  range  might  be 
achieved  by  maximizing  hypersonic  L/D, 

The  geometry  of  a  generalized  configuration  was  optimized  for  maximum  L/D  and  the 
following  characteristics  obtained: 

1,  Flat  bottom  lower  surface 

2.  Clipped  delta  planform  with  a  0.32  taper  ratio 

3.  Delta  dorsal  vertical  fins 

4,  (2:1)  elliptical  nose  and  leading  edges 

With  the  general  configurational  geometry  est-’hlished,  the  (L/D)  „  was  determined  for 

rndx 

the  composite  detailed  geometry  for  two  vehicle  weight  classes  (10,000  and  100,000 
pounds)  while  meeting  the  previously  established  11  constraints.  For  the  domain 
investigated,  these  numerical  results  show  that  low  aspect  ratios,  low  wing  loadings, 
and  large  scale  vehicles  produce  the  highest  L/D.  Also  the  temperature  constraint  for 
the  domain  investigated  is  the  most  severe  in  limiting  L/D.  Subsequent  to  the  temperature 
constraint,  the  effects  of  skin  friction  are  the  next  most  limiting.  Therefore,  with  the 
development  of  drastically  improved  materials  and  the  relaxation  of  the  thermal 
constraint,  higher  L/D  values  are  possible.  In  addition,  as  booster  payload  capability 

W 

increases  and  the  structural  designs  become  more  efficient  (reducing  -g-)  higher  values 
of  L/D  may  be  anticipated. 
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APPENDIX  I 

TRAJECTORY  ANALYSIS 


The  trajectory  must  be  computed  first  to  optimize  the  performance  of  a  vehicle. 

The  equations  of  motion  for  a  banked,  gilding,  re-entry  vehicle  are  listed  as  follows: 


mV  =  —  0  -  mg  sin  y 

(132) 

.  ,  m  V*  cos  r 

mV  y  -  L  cos  <p  4-  - ^ - - —  -  mg  cos  y 

(133) 

L  si  n  ^ 

mV 

COS  y 

(134) 

For  hypersonic  lifting  re-entry,  the  assumption  of  y  and  y  equal  zero  may  be  used.  Thus 

V  _  0 

g  W 

(135) 

L  V*  •  * 

—  cos  ^  s  1  -  - r-  ,  whtr#  Vg  s  g 

W  Vg 

"E 

(136) 

V'/'  L  .  , 

-  s  -  sin  ©  . 

g  w  ^ 

(137) 

The  range  equations  are  as  follows: 

Longiludinol  rong#  »  R*  =  / ''  V  cos 

Vi 

,  dV 
^  -7— 

V 

(138) 

0 

Loftral  rongc  =  Ry  ~  ~  J  V  tin 

Vi 

^  -r- 
V 

(139) 

where 

^  df  »  /  dV. 

Vi 

(140) 

When  the  equations  of  motion  (135)  and  (137)  are  used. 

'  V, 

(141) 

Integrating  Equation  (141)  for  constant  L/D  and  produces 

L  Vi 

«  —  sin  A  In  — “ 

^  D  ^  V 

(142) 
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Uaitiiv  liquations  (135).  (136).  nncl  (142)  In  Equations  (138)  and  (139)  results  In 

CO,  (t^i)  «ln  In  -5^)  d  (143) 


Vl 


=  /'  -F  co.^  sin  sin  in  d 


(144) 


Max'P’.iim  longitudinal  range  achieved  for  zero  bank  angle  (*  =  0)  in  Equation  (143) 
resuits  in  the  following; 


Vl 


/  \  .  Jz-  n  fi_  _5ii_  ^ 

'Re  'mo,  0  (  I- 2  0  ^  V^*  ^ 


(145) 


where  x  =  —  .  dummy  variable, 

c 

Figure  10  shows  a  plot  of  the  longitudinal  range  achieved  for  various  Initial  velocities. 

Vi  . 

The  integration  of  the  lateral  range  equation  is  more  Involved,  however.  For  1, 

c 


R  I 

(  rt)'  -q-  ^  /  rfr  *) 


(146) 


’E  '  “  o  I  -  K* 

Transformation  of  variables  permits  integration.  Let 

.  t 


and 


Hence 


y  =  —  In  * 


I  L  .  . 

o  »  —  -5-  sin  ^  . 


Y  -F  ^ 


9 

I 


•  (sin  0  y  ) 


4y 


r  _•  • 

s  J  s  ’  (  sin  0  y  )  £ 

0  k»  0 


,-ky 


dy 


(147) 
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In 


(^)  =  *  i 

t  n  =  I 

tlie  rL'^^ion  of  interest  (0  a  2),  tlte  series 


o*  +  n*  • 

may  l>e  approximated  as  follows; 


(149) 


i 

(^)Sin  2 

4> 

5.2 

1  +  0. 106 

(iy 

sin*^  1 

(150) 


Hy  means  of  differentiation,  the  optimum  Ixink  angle  may  be  obtained  to  achieve  maximum 
range  as  follows; 


0  '06  (t)*  •  (151) 

Figure  11  is  a  plot  of  optimum  bank  angle  as  a  function  of  L/D  as  obtained  from  Equation 

(151) .  ny  insertion  of  Equation  (151)  into  Equation  (150)  the  maximum  range  for  the 
optimum  bank  angle  is  obtained, 

(t)' 

5. 2 /it  0. 106  (i)' 

Figure  12  presents  maximum  lateral  range  versus  L/D  as  obtained  for  Equation  (152). 

It  is  apparent  from  Equation  (145)  or  Figure  10  that  maximum  longitudinal  range 
control  is  achieved  by  modulating  L/D  at  hypersonic  speeds.  As  apparent  from  Equation 

(152)  or  Figure  12,  tlie  maximum  lateral  range  control  can  be  achieved  by  maximizing 
L/D  during  re-entry.  Thus,  maximum  maneuverability  is  achieved  by  maximizing 
hypersonic  L/D. 
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APlMiNOIX  II 

III- AT  TRANSI'HIl  ANALYSIS 


ll-A.  Intrcxluction 

The  desipn  of  hyper  velocity  re-entry  vehicles  depends  vitally  upon  the  heat  transfer 
rates  and  temperatures  that  occur  at  critical  points  on  the  vehicle.  The  re-entry  vehicle 

2 

of  the  lifting  or  glide  type  experiences  much  less  severe  heating  rates  (tens  of  BTU/ft 

2 

sec.)  than  for  the  ballistic  type  (thousands  of  BTU/ft  sec.),  but  the  total  heat  transferred 
is  greater  because  of  the  extremely  long  flight  times  of  the  former  in  the  atmosphere 
(al)out  an  hour  as  opposed  to  a  few  seconds).  Because  of  this  fact  ablation  techniques 
have  not  proved  to  be  practical  and  a  radiation-cooled  structure  is  used.  In  lifting  re- 
Cii;  V  a  steady-state  thermal  condition  is  quickly  reached  in  which  the  influx  of  aero¬ 
dynamic  convective  heating  is  balanced  by  the  efflux  of  thermal  radiation  resulting  in 
an  equilibrium  temperature  of  a  surface.  This  temperature  is  maintained  below  the 
limit  of  the  surface  material  (approximately  3000®F)  by  carefully  designing  the  shape 
of  the  vehicle  and  designating  flight  paths  that  avoid  the  critical  heating  region. 

The  .stagnation  region  usually  experiences  the  highest  heating  rate  on  hypersonic 
vehicles.  Many  methods  are  available  for  predicting  laminar  heating  rates  at  the  stagna¬ 
tion  pojnt  and  these  will  be  Investigated  with  the  objective  of  ascertaining  the  simplest 
method"  that  possesses  sufficient  accuracy  for  the  heating  analysis. 

The  leading  edge  is  the  next  most  important  area  to  be  considered.  Transformation 
is  possible  to  convert  from  three  dimensional  to  two  dimensional  stagnation  heating. 

The  leading  edge  geometry  in  section  and  planform  dictates  the  heating  rate  distribution. 

Next,  the  analysis  of  the  lower  surface  aft  of  the  nose  and  leading  edge  must  be 
analyzed  for  both  laminar  and  turbulent  flow. 


II -B.  Stagnation  Point  Heating 

Several  of  the  available  methods  for  computing  stagnation  heat-transfer  rates  were 

w 

studied  and  computed  for  a  typical  lifting  re-entry  glide  path  —  =  100)  as  shown  in 
Figure  13,  The  methods  analyzed  are  listed  as  follows: 

1.  Sibulkln 

2.  Cohen  and  Reshotko 

3.  Lees,  approximate  and  exact 

4.  Uetra,  Kemp,  and  Riddell 

5.  Fay  and  Riddell 
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Figure  13.  Comparison  of  the  Stagnation-Point  Heating  Rate  for  Various  Methods 


The  simplest  method  (Lees  approximate)  also  Is  closest  to  the  numerical  average  of  all 
methods.  For  this  reason  Lees’  approximate  method  In  the  following  form  was  adopted: 


q  =  (15,5  )  2* 

s 

II-3.  Leading  Edge  Heating 


(153) 


The  unswept  leading  edge  may  be  obtained  simply  from  Lees’  stagnation  equation  for 
n  "  0.  The  conversion  for  sweep  as  shown  in  Reference  5  is  as  follows: 


where 


(  sin  A  cos  a  )  for  the  wing  landing  adga 


(154) 


(155a) 
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or 


Thus 


A.  =  Ar  +  <1  I  for  !*'•  Hn  laading  adg*. 


■^LE  =  ‘>S 


V  ZRle  } 


0.  s 


:ot  A. 


(ISSb) 


(156) 


11-4.  Lower  Surface  Heating 


The  nose  and  leading  edge  are  generally  laminar  at  critical  heating  for  “state  of  the 
art”  materials,  structures,  and  geometries.  Aft  of  the  leading  edge  local  transition  to 
turbulent  flow  may  occur  depending  upon  attitude,  altitude,  velocity,  surface  roughness, 
temperature,  and  geometry.  As  an  allowance  for  a  degree  of  conservatism,  the  lower 
surface  will  be  designed  for  either  laminar  or  turbulent  flow  at  critical  heating  depending 
on  which  is  higher. 


The  method  adopted  for  computing  heat  transfer  is  the  reference  enthalpy  technique 
first  expressed  by  Eckert  (Reference  6).  Schmidt  (References  7  and  8)  and  Hankey 
(Reference  5)  have  computed  laminar  and  turbulent  heat  transfer  to  flat  plates  for  re¬ 
entry  flight  conditions.  An  empirical  approximation  to  these  data  applicable  in  the 
critical  heating  regime  are  tabulated  as  follows- 


4 

Lam  FP 

4 

Turb  FP 


12.1 

.0,9 


/ — 5^* — / 1  _  JijL  )  (  p  *in*  a  009  a  ) 
\  (000  /  \  h,  / 


1H2.  (JU-I  U-^)  I/. 

*  0.2  '  1000  '  '  h  -  ' 


.  t  .0.9 

Sin  a  009  a  i 


(157) 

(158) 


11-5.  Critical  Heating 


For  equilibrium  glide  (Appendix  1)  in  which 


7 


w 

SCl 


(159) 


the  flight  conditions  for  which  maximum  heating  occurs  can  be  determined  since  q  is  pro- 


m  3 

portional  to  o  _  (where  m 
q  is  a  maximum  when 


■  0.8  for  turbulent  and  0.5  for  laminar  heating).  Therefore, 


dp  _3p 

dV.  ■  mV, 


(160) 


From  the  differentiation  of  the  equilibrium  glide  Equation  (159),  however, 

(,  _  JL.  )  ' 

IZ  '  V.  V,*  I 


(161) 
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When  liquations  (160)  and  (161)  are  combined,  the  velocity  at  critical  heating  may  be 
ascertained. 


^  lo. 


8165  lomlnor 
6831  turbuUnI 


or 


(  2  I ,  I  70  lam  I  nor 
^mox  (  17,710  turbulant 


(162) 


(163) 


Thus  the  critical  heating  rates  on  the  vehicle  during  equilibrium  glide  re-entry  may  be 
determined  by  substituting  Equations  (159)  and  (163)  into  Equations  (153),  (156),  (157), 


w 


aiv!  (158),  assuming  —  «  1.  Therefore , 
"s 


,  Ml  P-* 

(q,  )  =  8.025  (— - —  )  ; 

'  *  'mo*  '  R|g  SCl  ' 


laminar 


,  W  cos*Aa  V®  *  . 

(q,_)  =  5.675  ( - *-)  ,  lominor 


mqi 


Rle  SCl 


<\sJmo*‘  »  (tf) 

/  w  \*. 


lurbu  lanl 


(164) 

(165) 

(166) 

(167) 


For  thermal  equilibrium  conditions  in  which  the  convective  aerodynamic  heating  is 
balanced  by  the  surface  thermal  radiation  to  space,  the  temperature  may  be  computed 
in  the  following  manner; 


Thus 


'LE 


4 

s  €  O’ 

T* 

(168) 

3594  ( 

0.01  W 

1 

(169) 

SCl  Rn  * 

3297  ( 

o 

o 

< 

CO**  A«  \* 

(170' 

SCl 

R  ' 

'‘UE 

2(30 

,  0.1  w 

1 

^  . 

1  am  1  nar 

(171) 

^  X  S 
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**  =  2333  ^ — j - j  turbultnt  ^172) 

K  *  S 

liquations  (169)  and  (170)  may  be  used  to  determine  the  minimum  nose  and  leading  edge 
radii  for  specified  material  temperature  limits,  sweep  angle  and  glide  parameter 
W 

(?7r~).  Figure  14  is  a  plot  of  minimum  nose  radius  versus  temperature  limit  for  an 

f  =  0.8. 

liquations  (171)  and  (172)  may  be  used  to  determine  the  maximum  wing  loading  (W/S) 
for  a  given  material  temperature  limit.  Figure  15  is  a  plot  of  maximum  wing  loading 
versus  lower  surface  temperature  limit  for  an  t  of  0.8,  x  of  1  foot,  and  a  Cj^  =  2  sin2Q 

cos  o. 


11-6.  Material  Temperature  Limits 


High  temperature  materials  are  required  to  permit  successful  re-entry  of  glide 
vehicles.  The  group  of  materials  most  applicable  for  this  purpose  is  entitled  "refractory 
metals."  The  melting  points  of  some  of  the  most  applicable  refractory  materials  are 
tabulated  as  follows: 

Metal 

Melting  Point 

Tungsten 

6170*F 

Tantalum 

5425'F 

Molybdenum 

4739“F 

Columbium 

4380'’F 

Although  the  melting  temperatures  are  quite  high,  these  materials  all  suffer  drastically 
from  unsatisfactory  oxidation  resistance  at  temperatures  in  excess  of  2{XX)*F.  The  most 
desirable  method  for  overcoming  the  disadvantages  of  poor  oxidation  resistance  is 
through  alloying.  Since  little  improvement  in  the  oxidation  resistance  of  alloys  has  been 
accomplished  in  the  past  years,  it  is  generally  conceded  chat  coatings  are  required  to 
protect  structural  refractory-metal  alloys  from  oxidation.  In  addition,  the  bare  metals 
possess  low  values  of  emlttance;  hence,  a  material  coating  with  high  emittance  is  desired 
for  application  to  radiation-cooled  glide  re-entry  vcliicles.  The  problem  is  reduced  to  one 
of  finding  a  high  emittance  coating  that  can  withstand  a  wide  variation  in  temperature 
and  pressure  in  an  oxidizing  atmosphere  for  prolonged  periods.  Reference  9  indicates 
that  coatings  have  been  developed  with  many  hours  of  life  at  2S00*F  and  possess  a  cap¬ 
ability  at  3000*F  for  short  durations. 

Since  the  nose  of  the  vehicle  represents  a  small  portion  of  the  weight  yet  can  contrilHite 
a  large  amount  of  drag  and  appreciably  affect  the  maximum  L/0,  a  nose  cap  should  be 
used  with  a  greater  temperature  capability  than  the  remainder  of  the  vehicle.  The  coated 
carbides,  graphites,  and  various  oxides  possess  high  allowable  thermal  limits  (with 
reduced  strength  over  refractory  metals)  and  represent  potential  nose  cap  material.  The 
present  "state  of  the  art”  thermal  capability  of  nose  cap  material  is  in  the  vicinity  of 4000*F. 
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Figure  14.  MlnimuiK  lUdlu*  for  Allowable  Noee  Temperature 
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APPENDIX  III 

HYPERSONIC  AERODYNAMIC  CHARACTERISTICS 


Ill-A.  Introtluction 

There  is  not  available  presently  an  all  inclusive  compilation  of  procedures  for 
computing  re-entry  vehicle  force  and  moment  coefficients  at  hypersonic  speeds.  Simple 
techniques  are  developed  in  this  appendix  for  rapidly  determining  the  aerodynamic 
pressures  and  resultant  six-component  force  and  moment  characteristics  for  hyper¬ 
sonic  glide  vehicles.  The  methods  are  applicable  between  10*  and  30*  angle  of  attack 
since  this  covers  the  (L/D)  range  of  Interest  (from  4  to  1).  For  extremely  low  or 

exiretnely  high  angles  of  attack,  prediction  of  the  aerodynamic  characteristics  becomes 
mo'  w  complex  and  need  not  be  investigated  here.  Techniques  herein  are  further  limited 
to  the  continuum  flow  regime  in  that  the  aerodynamic  characteristics  are  only  of  primary 
interest  (for  stability  and  performance)  when  the  dynamic  pressure  Is  significant.  At 
low  dynamic  pressures  (q  10  pounds  per  square  foot)  reaction  controls  are  required 
for  stability  and  control  of  space  vehicles.  At  orbital  speed  a  dynamic  pressure  of  10 
pounds  per  square  foot  occurs  at  an  altitude  of  about  265,000  feet  for  which  the  mean 
free  path  is  about  0.01  foot.  The  Knudsen  number  based  upon  a  10-foot  chord  would  be 

.*1 

10  indicating  that  a  continuum  flow  model  may  be  used. 

Expressions  for  the  pressure  distribution  over  simple  shapes  (spheres,  cylinders,  and 
planar  surfaces)  are  first  derived.  A  generalized  configuration  is  defined  composed  of 
the  preceding  simple  elements,  and  an  analytical  expression  for  the  aerodynamic  co-. 
efficients  resulting  from  integration  of  the  surface  pressures  is  presented.  Finally,  a 
summary  of  the  six-component  aerodynamic  characteristics  is  tabulated  convenient  for 
programming  on  high  speed  computers. 


II1-F3.  Pressure  Distribution  Theory 

For  rapid  prediction  of  aerodynamic  forces  and  pressures  on  hypersonic  aerodynamic 
vehicles,  simple,  basic  theory  Is  desirable.  Perhaps  the  most  popular  theory  for  this 
type  of  analysis  is  Newtonian  Impact  Theory  and  Its  various  modifications  (References 
10  through  19). 

The  Newtonian  flow  concept  assumes  that  the  free-stream  gas,  upon  striking  a  surface, 
loses  its  component  of  momentum  normal  to  that  sifrface  and  then  moves  along  the 
surface  with  Its  tangential  component  of  momentum  unchanged  (Inelastic  collision).  From 
this  assumption,  the  pressure  coefficient  Is  defined  by: 

Cp  »  2  tin's  073) 

where  fl  -  local  flow  inclination  with  the  free-stream  (angle  measured  between  the 
velocity  vector  and  the  surface)  and 
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■■  relarions  are  modified  for  the  hypersonic 
ICS  the  Ixxly  slope  and  the  identical  result, 
basic  validity  of  the  Newtonian  Flow  Theory 


Lees,  in  Refen  ■  i  ;  t  21  :\1  Ncwfonian  Theory  as: 


where 


5in *  8 


(174) 


1  i  '  pressure  coefficient. 


A  further  mtxlific. It i'  1  i-i*  ;  •  C‘ 
Nc'vronian  Theory 

(' 

1  , 

'  r  5  I 


Reference  16  is  called  the  Generalized 


sin*  8„ 


(175) 


This  form  of  tiie  Newa  n  u  1  1  ;c  •  asetni  for  pointed  shapes  such  as  tangent  ogives. 


In  both  of  rlie  preceding,  '•  x  i(  :  a  :ons,  ’he  basic  form  of  the  equation 

=  X  sin*  8  (176) 

is  retained.  The  k -factor  refleci  s  :ae  change  in  pressure  coefficient  because  of  such 
factors  as  initial  nose  slope,  Maci  number,  and  change  in  gas  composition,  while  the 

8in'^6  term  accounts  for  the  surface  orientation  and  geometry.  Newtonian  Flow  Theory 
in  its  basic  form,  l-quation  (P3),  us  applicable  only  to  hypersonic  flow  over  highly  inclined 
surfaces  (Reference  Id);  however,  extension  of  the  theory  is  possible  if  the  form  in 
Equation  (176)  is  retained  and  the  k-factor  adjusted.  By  the  use  of  the  theories  of  known 
applicability  and  experimental  data,  the  k-factor  for  simple  geometries  is  determined 
herein. 

The  chief  advantage  in  retaining  the  expression  of  the  pressure  coefficient  in  the  form 
of  Equation  (176)  is  found  in  the  determination  of  force  coefficients  in  that  the  required 
integrations  are  made  relatively  simple.  A  further  advantage  exists  whereby  the  force 
coefficient  is  directly  proportional  to  the  k-factor.  Therefore,  the  variation  of  k-factor 
over  a  range  of  flight  conditions  directly  defines  the  variation  of  the  force  coefficients. 


Ill-D-I.  Blunt  Surfaces  In  Impact  Flow 

Variation  of  the  k-factor  for  the  stagnation  region  of  a  surface  is  illustrated  in  Figure 
16.  In  this  figure,  the  variation  of  k  for  a  real  gas  as  computed  by  Feldman,  Reference 
20,  is  shown  for  a  range  of  Mach  numbers  and  altitudes.  Also  shown  is  the  variation  of  k 
for  an  ideal  gas,  as  defined  In  Reference  21  and  determined  from  the  following  relation¬ 
ship; 
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II  I  i  n  M;ii.li  nunilx;r.  Figure  16  can  be  used  to  determine  the 
!  ic'i  ■  .  Note  that  real  gas  k-factor  values  are  seen  to  be  rather 

I  "ar  i:i-  chosen. 


I'lie  api'Mc.il  1  ’  (  f  Vim  i'  ed  \i'v/tonian  Theory  on  hemispheres  and  hemicylinders  is 
ilUisirated  in  !  carj  I  ’hi  figure  illustrates  that  Fquation  (176),  where  k  is  determined 
frot  1  I'igure  !  .  i:  i  i a,  evaluating  pressures  and  forces  on  these  suffaces. 


Figure  16,  Variation  of  Stagnation  Pressure  Coefficient  with  Match  Number 


In  the  analyzation  of  swept  leading  edges,  cross-flow  theory  Is  used.  This  reduces 
the  complex  three-dimensional  problem  to  a  simpler  one  In  two  dimensions.  Using  cross- 
flow  theory  Equation  (176)  becomes 

Cp^  «  k  co»*A,  (178) 

where  is  the  effective  sweep  angle  of  the  leading  edge,  see  Equation  (196). 
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Figure  17.  Applicability  of  Newtonian  Theory  to  Hemispheres  and  Hemicylinders 


ni-B-2.  Flat  Surfaces  In  Impact  Flow 

Where  the  Modified  Newtonian  Theory  Is  significantly  in  error,  and  the  surface  under 
consideration  is  large,  some  correction  must  be  applied.  An  example  of  such  a  situation 
is  the  large  lower  surface  of  a  blunt  swept  delta  wing.  Figure  18  shows  the  change  in 
pressure  coefficients  on  a  blunt  wedge  due  to  change  in  angle  of  attack  (Reference  22). 
For  angles  of  attack  between  5*  and  15®,  the  "tangent  wedge”  theory  appears  representa¬ 
tive  of  the  mean  of  the  data.  For  angles  of  attack  above  15®,  the  flow  appears  to  change  in 
nature  and  approaches  Newtonian  values  until  at  90®  the  flow  stagnates  at  C  . 
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As  no  sinj’le  simple  tlioory  wouKI  predict  the  change  in  the  nature  of  the  flow  across 
tlie  angle-of-attack  spectrum,  several  empirical  correlation  schemes  were  attempted. 
The  most  successful  of  tlic  schemes  considered  is  the  one  shown  in  Figure  19.  The 
faired  straight  line  through  the  data  has  the  equation 

k  =  1. 9P  +  0.  21  col  8  .  (179) 


8  (d«gr«M) 

Figure  18.  Low«r-6urfsoe  Preiaure  Coettlolents 
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Figure  19.  Empirical  Correlation  of  Lower-Surface  Presaure  Coefficients 


For  surface  inclinations  below  10°,  the  representation  in  Figures  18  and  19  rapidly 
loses  significance  because  at  6  =  0  for  any  Mach  number,  k  approaches  » ,  which  requires 
a  change  in  the  expression  for  Cp.  Interaction  and  induced  pressure  effects  also  become 

dominant  at  low  angles  of  attack,  which  require  a  change  In  the  analytical  procedure.  For 
these  reasons,  values  of  at  much  below  10*  will  not  be  considered. 

More  comprehensive  analyses  must  be  applied  only  when  the  vehicle  in  question  is 
considered  to  fly  at  low  angles  of  attack.  As  presently  anticipated,  only  gliders  with  L/D 
greater  than  4  will  fly  at  angles  of  attack  less  than  10*.  These  vehicles  will  require 
refined  analysis  at  low  angles  of  attack. 
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llI-n-3.  Iruluct-'il  Pressures 

.1  blunted  wed^e  at  angle  of  attack  produces  a  curved  shock  wave  near  the  leading  edge. 
Farther  away  from  the  leading  edge,  the  shock  wave  isymptoclcally  approaches  that  of  a 
sharp  wedge.  Hecause  of  this  shock  curvature,  induced  pressures  commence  at  the 
leading  edge  shoulder  and  diminish  as  the  shock  curvature  reduces  (References  23  and  24). 
Techniques  are  available  (e.g.  Relotserkowski,  Vaglio-Laurin,  etc.)  that  provide  numerical 
solutions  of  the  flow  field  from  which  pressure  distributions  may  be  ascertained.  However, 
since  numerical  results  are  not  yet  available  in  quantity,  empirical  correlation  of  the 
induced  pressures  was  attempted. 

From  experimental  pressure  distributions,  we  determined  that  a  linear  variation  in 
AC  from  /.ero  at  the  centerline  to  AC  at  the  shoulder  would  adequately  account  for  In- 
^  ^sh 

duced  pressure  distribution  effects  on  delta  wings  with  relatively  small  leading-edge  radii. 
Frnn  ■r.e  leading-edge  pressure  distribution  of  Reference  22,  we  observed  that  the  shoulder 
pressure  may  he  ascertained  by  shifting  the  flow  angle  for  the  Newtonian  pressure  coefficient 
by  a  small  value,  {  ,  so  that 

Cp  =  k  »in*  (a  +  <  1  .  (180) 

The  induced  pressure  coefficients  at  the  leading  edge  shoulder  are  presented  in  Figure 
20  for  various  a  values.  From  this  figure,  it  Is  apparent  that  the  form  of  Equation  (180) 
adequately  accounts  for  Induced  pressure  if  a  value  of  (  of  14*  is  included  while  the  un¬ 
modified  Newtonian  (f  =  0)  is  quite  inadequate.  For  swept  leading  edges 

Cp  ^  =  h  cot*  A,  tin*  (Op  +«  1  (181) 

where  and  are  defined  by  Equations  (202)  and  (203).  When  the  identity  is  used, 

cot*  Ap  tin*  Op  =  tin*  a  cos*  0  .  (182) 

Equation  (181)  may  be  expanded  for  small  values  of  jS  and  f  as  follows; 

Cp  =  k  tin*  a  (  I  +  Zt  cot  Op  )  .  (183) 

The  difference  in  shoulder  pressure  values  is  required  in  the  determination  of  the 
rolling  moment  due  to  yaw  of  delta  wings.  This  pressure  difference  can  be  obtained  by 
using  Equation  (183)  for  the  right  and  left  shoulders  as  follows: 

A  Cp  e  2  k  €  tin*  a  ( cot  Op  -  col  Op  ) .  (184) 

R-L  ***  R  L 

For  the  geometry  of  interest  ir  this  study,  we  determined  from  Equation  (203)  that  ' 


A  col  Op  •»  4  /9  . 
R-L 
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Figure  20.  Induced  Pressure  Correlation  of  Blunt  Leading  Edges 


Therefore, 


A  Cp..  »  Xi  /3  *in*  a 

R-L 
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Ill-H-'l.  Snfl'.K'c  S)-- ill  I'rirtioii 


Several  lUL'ilioJs  have  iK'cn  prost  iiteil  in  the  literature  for  the  prediction  of  laminar 
skin  friction  on  flai  plates  in  8,iper:.onie  and  hypersonic  flow.  These  methods  are 
usually  ci'inplo;  anti  lalwiriotis  to  apply  and  often  require  evaluation  of  the  flow  field 
iK'hind  the  compression  shock  wave.  One  of  the  more  successful  methods  is  the  Reference 
linthalpy  NUnhod  as  tlefined  in  Reierenee  6.  .I.F.  Schmidt,  in  Reference  7,  has  used  thi.s 
methtxi  to  predict  skin  friction  over  a  flat  plate  for  a  wide  range  of  flight  conditions. 

Since  the  data  arv  ^iven  in  tertiis  of  flight  ptirameters  instead  of  flow  parameters  at  the 
edge  of  the  Ivnindary  layer,  this  reference  i.s  ideal.  For  further  simplification  of  the 
pre  '  ;  am  of  laminar  skin  frieiion,  an  empirical  curve  was  fitted  to  the  data  presented 
by  Sclanitlt.  This  single  line  is  represented  hy  the  equation 


lam 


yw: 


=046  cos  a  +4  65 


10  ,  OOC 


sin  a 


z.z 

cos  a 


(187) 


where  o  is  the  surface  angle  of  attack  and 


C 


I 


:  locol  shsor  stress. 


(188) 


riiis  ei|ua',ioii  deviates  no  more  than  20  percent  from  the  data  presented  by  Schmidt  for 
low-altitude  high-angle-of-attack  flight  and  i.s  m.orc  accurate  (approximately  by  ±10 
percent)  for  the  rest  of  the  altitude  and  angle-of-attack  range.  This  accuracy  is 
sufficient  for  preliminary  design  as  long  as  the  vehicle  is  in  continuum  flow  where 
shear  stresses  do  not  dominate,  A  comparison  of  the  preceding  equation  with  the  method 
of  Schmidt  is  shown  in  Figure  21.  This  equation  is  seen  to  be  applicable  for  the  majority 
of  the  re-entry  flight  corridor. 


As  is  the  case  with  laminar  skin  friction, many  techniques  are  defined  in  the  open 
literature  for  the  prediction  of  turbulent  skin  friction  on  flat  plates  in  supersonic  and 
hypersonic  flows.  Unfortunately,  these  methods  are  even  more  complex  than  those  for 
laminar  flow  and  rely  cn  emipirical  rclation.s  derived  from  wind-tunnel  or  flight-test 
data.  Analysis  is  often  further  complicated  by  the  need  to  evaluate  the  flow  field 
properties  outside  the  boundary  layer.  Schmidt,  in  Reference  8,  has  used  the  reference- 
enthalpy  technique  to  construct  a  comprehensive  series  of  graphs  for  the  prediction  of 
turbulent  .skin  friction  on  a  flat  plate  for  a  wide  range  of  flight  parameters.  As  in  the 
laminar  study  hy  Schmidt,  the  data  are  given  in  terms  of  flight  parameters,  not  local 
flow  parameters,  and  are,  therefore,  quite  useful  in  preliminary  design.  An  empirical 
curve  fit  of  Schmidt’s  data  resulted  in  the  following  equation: 

C,  (R.  1°*  =  0.  Os  8  (sin  4.5  a  1+0.70  - ^ —  COS**®  O  sin''*a 

'lufb  "«  10,000 


The  degree  lo  which  this  equation  covers  the  analysis  of  Schmidt  is  seen  in  Figure  22. 
rtilH  close  comparison  proves  the  equation  sufficiently  accurate  for  preliminary  design 
purposes. 
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APfROX 


Figure  22,  The  Compariion  o(  the  Empirically  Oetermiaed  Turbulent-Skin  Frlotlon  Equttlou 


III-C,  Determination  of  Aerodynamic  Characteristics 

The  purpose  of  this  section  Is  to  apply  the  previously  developed  theoretical  and 
empirical  relations  to  the  various  basic  components  of  a  generalized  hypersonic  llftini^ 
re-entry  vnlilcle  and,  considering  the  relative  location  of  the  components,  describe  the 
aerodynamic  characteristics  of  the  configuration,  Dependence  of  the  aerodynamic 
characteristics  on  surface  geometry  and  orientation  is  reflected  in  the  Newtonian 
format 

Cj  ■  h  lia*#  (176) 

as  developed  In  the  previous  section.  Determination  of  6  and  Integration  of  for  the 
required  aerodynamic  characteriatica  of  surfaces  with  geometries  not  easily  expressed 
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analytically  often  can  be  done  precisely  only  by  numerical  methods.  This  complexity  is 
certainly  to  be  avoided  whenever  possible.  For  this  reason  only  simple  geometric  shapes 
will  lie  analyzed.  Quite  accurate  analysis  of  complex  geometric  shapes  can  be  made  by 
proper  selection  and  combination  of  the  simple  shapes  covered  herein. 

As  important  as  the  definition  of  the  components  is  their  location  relative  to  the  center 
of  gravity  of  the  vehicle  in  question.  Therefore,  the  geon.etry  and  the  coordinate  system 
of  the  generalized  configuration  to  be  analyzed  must  be  established.  Efficient,  lifting 
hypersonic  vehicles  possess  a  flat  bottom  surface  with  a  highly  swept  planform  to  maximize 
Cj  aivl  L/D.  The  nose  and  leading  edges  must  be  blunt  (i.e.,  circular,  elliptical,  etc.)  to 

withntand  the  heating  environment.  Vertical  fine  are  required  to  provide  directional  sta¬ 
bility  and  are  most  practical  and  effective  outboard  and  in  the  most  rearward  position. 

The  body  is  located  on  the  upfier  surface  unexposed  to  the  high  heating.  Figure  23  is  a 
three-view  sketch  of  the  generalized  configuration  to  be  considered  showing  the  general 
dimensions  of  the  overall  configuration,  the  various  components,  and  the  location  of  the 
ft  -ter  of  gravity. 


rifuiw  2).  Ttw  OaiwraUaed  Conf ttpuratioa  Oaomatry 
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With  a  j;onerali;;ed  configuration  defineil,  the  aerodynamic  charucterlstlce  can  be 
analyzed  in  terniB  of  the  general  dimensions  of  the  configuration.  Defining  these  general 
dimensions  for  a  specific  configuration  permits  rapid  determination  of  the  aerodynamic 
characteristics  of  that  configuration. 

1'ho  aerodynamic  coefficients  derived  herein,  and  the  associated  sign  convention,  are 
defined  in  tlie  list  of  symbols.  The  symbol  and  sign  convention  used  by  the  NASA  has 
been  used  unless  otherwise  noted. 

IXirlvation  of  tiie  coefficients  for  the  different  components,  in  genera!,  followed  three 
basic  steps: 

1.  Determination  of  a  simple  representative  shape 

2.  Integration  of  the  theoretical  pressure  coefficients  over  that  simple  shape  and 
evaluation  of  the  basic  force  coefficients 

2.  I>jtermi nation  of  appropriate  moment  coefficients  through  consideration  of  the 
displacement  distances  of  the  component  from  the  vehicle  center  of  gravity. 

Six  component  characteristics  were  defined  for  each  configuration  component.  The 
lengthy  details  of  the  integrations  required  to  derive  the  aerodynamic  force  coefficients 
were  omitted  and  only  the  results  presented. 


III-D.  Summary  of  [>;rived  Components 
HI  D-1.  General 

A  summary  of  all  the  derived  coefficients  are  tabulated  for  convenient  reference.  Only 
the  basic  equations  arc  given.  The  terms  such  as  k-factor,  a,  A^,  a^,  etc.  in  the 

equations  must  be  made  to  apply  to  the  configuration  component  that  is  evaluated. 

The  lateral  and  directional  coefficients  were  determined  to  be  linear  with  fi  for  small 
angles  and  hence  only  the  derivatives  are  presented.  Numerical  values  of  the  derivatives 
are  In  radian  measure. 
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in-n-2.  Nohc  liquations 


w  Hn 

sin  a  ( 1  +  cos  a  ) 

(190) 

'-N 

AS 

w  Rfj*  k|^ 

( 1  4  cos  a  )* 

(191) 

AS 

2 

Cm 

-  C  - 

Ca 

c 

(192) 

Ca 

(193) 

cos  a 

C4 

(194) 

cos  a  b 

/« 

1 

0 

> 

*N 

(195) 

cos  a 

b 

All  derivatives  and  all  angles  are  in  radians. 
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ni-l)-3.  Lcjuling  lidgc  liquations 

These  equations  are  for  a  pair  of  leading  edge  elements. 


t  “^LE  \  . 

\  35  )  ''le  *•'>  O  (co»A,+cosA  cos  a  ) 


^  cos  A  (cos  A^  +  cosA  cos  a  )* 


‘LE 


=  Cn  -Ca 

c  c 


LE 


'  'i 


/  4R|  c  •^IF  \  *  t 

*  ~  V — - /  •'i.E 


tiLE  •'"A  (1+ cos  a,)  [  — ^  sin  A  cos  A*  (14-cos  a,) — ^ 


sin  A  cos  A,(I4-cos  a,)*( 


* /•*LE 


sin 


A+  cos  a) 


(sin  Af)^  =  sin  A  cos  a  cos  0  7  cos  A  sin  /3 
L 

«  .  .  A  ^  sinA  a 

(cot  a,)_  =  cos  A  col  a  ±  -  lonfl 

•  R  sin  a 

L 

for  one  leading  edge 


(196) 

(197) 

(198) 

(199) 

sin  a] 

(200) 
(201) 

(202) 

(203) 
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Ul-I)-4.  Lower  Surface  L(|uutions 


=  kjj  (-^)  ‘la  a 


(204) 


C*  =  <3 


S  0.45  cos  a  +4.65 


10,000 


sina  cos*‘*a 


(laminar)  (205) 


,  S  .  0.048  sin  (4.5  a) +  0.70 


(“s-) 


{^s 


10.000 

t  f 


*.*S  .  I.5_ 

cos  a  sin  a 


(turbulant)  (206) 


(207) 


0 


Ca 

cos  a 


(208) 


^A  JiL 
cos  q  b 


kj  sin^a 


IlLE 

9b 


k|  s  3.8 


C 


"0 


»L 

COSO  b 


where 


I  ♦  n 

C  -  [-^=^ 

n  (  I  +  n  )  '•  I  -  m 

_  1  0 . 5  lominor 
"  '  I  0.8  (urbulont 

m  :  plonform  topsr  ratio. 


(209) 

(210) 
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in-D-5.  Vertical  Fin  FquatlonB* 


Cn 

8Rp-tp  kLE  t  .  .  .  .  . 

“  “  •  2s  (Ap-fa)  »in  Ap 

(211) 

Ca 

Sp  -3  2,  ®^F^F  ^LE  2  a  a 

=  2kc  —  IX  co«  a)+ - cos‘(Ar  +  a)  co«  Ap 

(212) 

Cfn 

-  -  2kp  (X®  cos*  a )  - cos*I  Ap  +  a  ) 

sin  Ap  +  COS  Ap  ] 

(213) 

Sp 

=  -4kF  (X  cos  O  ) 

(214) 

Sc  /  *F  V 

*  -4kp  -|-  (Xcosa  )  (-  -g-) 

(215) 

,  4kp  (Xcos  a)  (- +  xJ^) 

(216) 

Sp 

for  on«  fin 

-^F 

for  on*  fin 

•The  equations  are  for  a  pair  of  fins. 
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UI-l)-6.  lilriuations  for  Transfer  of  Moment  Re^'erence  Center 


Ax  Az 

ACm  -  ~  — 

C  r. 


Az 


AC 


b 

Ax 


(217) 

(218) 

(219) 


in-l)-7.  Transformation  liquations  From  Body  to  Wind  Axes 

Cl  *  Cm  cos  a  -C^  sm  a 
Co  *  Cm  sin  ot  +CaCos  a 


(220) 

(221) 
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AI’Pi:NniX  IV 
STAniLlTY  ANALYSIS 


The  dynamics  of  a  vehicle  must  be  such  that  the  handling  qualities  are  within  the 
servo-response  capabilities  of  the  pilot.  The  short  period  modes  are  those  of  primary 
concern.  The  dynamical  equations  of  motion  may  be  obtained  by  analyzing  the  moment 
about  the  oscillating  veliicle  as  follows; 


where 


T  = 

d  14 

(222' 

dl 

H 

=  3 

;  •  <£ 

(223) 

H 

K 

o 

0 

I. 

= 

0  ly 

0 

(  pr  Incipal  a««s  ) 

f224) 

0  0 

I* 

W 

-  o± 

+  i^k_ 

(225) 

dH 

8  H 

-  +  u> 

X  H 

(226) 

dl 

8t 

and  where  -A  denotes  differentiation  of  only  the  scalar  terms.  Thus,  when  Equations  (222) 
ot 

through  (226)  are  combined, 

I  =  [  +  (I,  -ly)  0  ^  ]  1 

+  [l,  e  4-  I  I,  -  ]  i 

+[i,  f  .+  <  ly  -  I*  >  ^  ^  ]  1  • 

For  a  pure  aerodynamic  torque  with  no  damping 

T  =  (qSb)C^i+(qSc)  +  (  qS  b  )  . 

The  preceding  static  aerodynamic  moments  may  be  linearized  as  follows; 

Ci  ’  ^Jlp 


(227) 

(228) 


Cflfi  »  Cm  Aa 


Cn  »  CbQ  . 


(229) 

(230) 

(231) 
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Ihe  aoroilymunic  damping  of  vehicles  at  hypersonic  speeds  is  virtually  nonexistent  in 

that  the  reduced  frequency  (•^^)  is  small.  For  example,  a  wing  with  a  50-foot  root  chord 

oscillating  at  1  cycle  per  second  at  Mach  Number  20  would  produce  a  reduced  frequency 

of  alx)ut  10  Negligible  damping  occurs  for  such  a  condition  and  the  dynamic  stability 
derivatives  may  be  omitted.  It  now  remains  to  relate  a  and  fi  to  6,  i(>,  and 

’I’he  velocity  vector  defined  in  terms  of  the  liuler  angles  for  the  pitch-yaw-roll 
rotational  sequence  (0  -  i/i  -  0)  is  as  follows  (Reference  25); 

y 

=  cos  cosdj^-f(sin(^sind  -cos^sini/^cosdl) 

-t-(cos^sind  +  cos  0  tin}//  tini/>)il^.  (232) 

Fr“  small  angles  of  1/1  and  which  occur  as  perturbations  from  a  steady  level  flight 
condition. 


cos  9  i  +  i  <f)  sin  6  -  \/i  cos  0  )  +  sin  0  k_. 


f233) 


The  angle  of  attack  is  defined  as  follows; 

*  .X  :  L.  »  ^ 


ton  a  s 


cos 


0 


=  top  0. 


Thus 


and 


a  :  0 


(234) 


(235) 


a  =  0  (236) 

The  angle  of  sideslip,  |J,  measured  in  the  wind  reference  plane  is  defined  as  follows; 


sin  /3  - 


s  ~  •  j  s  ^  tin  9  -  }//  cot  9 


|V|  1Y|  ^ 

for  small  sideslip  angles,  and  since  6  =  a  (Equation  235), 

(3  =  <f>  tin  a  -  }//  cos  a 


(237) 


(238) 


and 

13  -  ^  sin  o  -  cos  a  tor  fixed  a .  (239) 

The  longitudinal  short -period  mode  for  a  hypersonic  vehicle  may  be  obtained  by 
setting  ^  and  ii  equal  to  zero  In  Equations  (227),  (228),  and  (230)  as  follows: 


ly  0  »  (d  S  c  )  A  a  ;  (-40) 

'  01 
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l«)wcvt>r,  i)  =  ft  from  llfiuation  (236);  therefore 

/  <1  S  F  \ 

a  (— -  )  A  a  =  0  ,  (211) 

ly 

wliich  is  the  ef|uation  for  simple  harmonic  motion  in  a  with  angular  frequency 


q  S  c 


(242) 


The  lateral  oscillation  short  period  mode  (Dutch  Roll)  may  be  obtained  by  setting  6 
equal  to  zero  in  liquations  (227),  (228).  (229),  and  (230).  See  Reference  26. 

I'herefore, 


^  =  (  q  Sb  CJp  )  A  fi 

Ij  (jit  =  I  q  Sb  Cnfj  )  A  0. 

Using  l:<iuations  (243)  and  (244)  to  eliminate  i  and  il  from  Equation  (2.39)  produces 

0  s  ^  ■  sin  a  -  — ^  cos  a  ^  ( q  Sb)  A /Q  , 

which  is  the  equation  for  simple  harmonic  motion  in  with  angular  frequency 

2  ,  Cnfl  . 

ton  =  q  Sb  (  — -  cos  a  ~  sin  a  j. 

P  '  I 2  Ik 


(243) 

(244) 


(245) 


(246) 


Pilot-rated  flight-simulator  studies  have  indicated  that  the  most  desirable  operating 
conditions  are  when  the  short  period  mode  of  the  vehicle  has  the  following  characteristics; 


frequency,  f  0.7 


f  >  0.4 
\  <-  1.0 


cycles  per  second 


damping  ratio,  C  -  0.7 


>  0.4 
<  1.0 


Oscillations  of  this  type  are  within  man’s  servo-response  characteristics  yet  the 
vehicle  possesses  satisfactory  sensitivity  to  control  inputs  (neither  too  sluggish  nor  too 
sensitive).  Present  day  aircraft  attain  these  handling  qualities  cither  by  aerodynamic 
means  or  artificially  by  "adaptive  control"  features.  Hypersoiiically,  negligible  aero¬ 
dynamic  damping  occurs  (C  —  0  for  the  emergency  situation  of  artificial  augmentation 
inoperative),  which  ncccHSltatcH  changing  the  handling  quality  criteria.  Since  zero  damping 
Implies  pcricKilc  motion  with  consiunt  umplitudc,  the  pilot  must  damp  out  oscillations 
by  "out-of-phase"  control  modulation.  For  accomplishment  of  this,  the  oscillation  must 
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Ik'  i»f  a  lonp  pet  UkI  (P  ■>  10  Koconds)  to  permit  sufficient  time  for  pilot  phase  control. 
I’lierefore,  with  zero  damping,  near  zero  frequency  is  required.  This  condition  is  not 
comfortahlo  to  fly  hut  considered  acceptable  through  flight  simulator  training. 

I'or  and  etiuai  to  zero,  then  from  iiquations  (242)  and  (243) 


Cffia  =  0 

(247> 

C„^  =  0 

(248) 

=  0. 

(249) 

lienee,  tiie  preceding  tiiree  conditions  are  the  hypersonic  stability  constraints  for 
flight  at  maximum  \  ./\). 
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AIMMiNDIX  V 
LOW  SLIilit)  ANAL  YSIS 


The  low-speed  flying  requirements  impose  several  constraiuis  on  a  hypersonic 
lifting  vehicle.  First,  the  vehicle  must  be  able  to  execute  a  satisfactory  horizontal 
landing  and,  secondly,  it  must  possess  acceptable  stability.  Low-speed  performance 
requirements  are  not  considered  important  since  most  of  the  maneuverability  Is  accom¬ 
plished  at  hyi^ersonic  speeds,  although  some  low-speed  maneuverability  is  a  fall-out 
from  the  landing  constraint.  The  point  is  that  the  vehicle  will  not  be  designed  for  good 
low -speed  performance  but  must  accept  the  amount  available. 


V-A.  L.anding 


The  method  presently  deemed  most  desirable  for  landing  hypersonic  aircraft  is  the 
“aiming  point"  technique  used  in  the  X-15.  In  this  method  the  pilot  dives  the  vehicle  at 
a  point  on  the  earth’s  surface  several  thousand  feet  before  the  runway,  then  executes 
a  flare  at  some  predetermined  altitude  to  a  shallow  glide  angle,  and  decelerates  by 
increasing  the  angle  of  attack  until  touchdown  at  some  preselected  speed.  By  the  use 
of  this  “Dynamic  Approach"  technique,  vehicles  of  low  subsonic  L/D  may  be  safely 
landed.  In  this  maneuver  the  highest  lift  coefficient  is  required  at  touchdown  (Reference 
27).  For  y  =  0 


(C.l 


t  'ma* 


ill 

'>T0 


(250) 


The  lift  coefficient  for  low  aspect  ratio  (/P)  airfoils  may  be  estimated  by  the  theory 
of  Jones  (Reference  28) 


C 


L 


a . 


(251) 


When  the  last  two  equations  are  combined. 


W  .  jr 

S  *2 


'TD 


(252) 


For  touchdown  velocities  less  than  about  200  knots  and  for  angles  of  attack  less  than  15° 


S  ^ 


i  50  pounds  p*r  squoro  foot. 


(253) 


Hence  the  preceding  relationship  between  wing  loading  and  aspect  ratio  becomes  the 
landing  constraint. 
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V-n.  Low  Speed  Stal)ility 

The  low-speed  directional  stability  constraint  Is  used  to  determine  the  size  of  the 
vertical  fin.  The  yawing  moment  for  two  dorsal  fins  may  be  expressed  as  follows; 


or  in  coefficient  form 


2Sf  \/“*F 


For  low  aspect  ratio  fin  planform  (see  Equation  251), 

•JCl  tt 


da 


F  ' 


(254) 


(255) 


(256) 


For  the  center  of  gravity  near  the  centroid  of  the  planform  (0.65  c)  and  triangular 
type  fin  planforms 


Also  for  delta  planforms 


and 


-  (0.35  c  1  . 


C  ^  A  A 

-  .  —  A 


3  2  cot  Af  . 

Thus  by  substitution  of  Equations  (256)  through  (259)  into  Equation  (255) 

O-  /  Sf  \  A 
■  3  \  S  /  ton  Ap  * 

For  clipped  deltas  of  taper  ratio,  m, 

m*  cot  Ap 


Sf 

s 


2  (I-  m*  )  cot  A 


Therefore  Equation  (260)  becomes 

'  3  V  m*  S  ^  • 


(257) 

(258) 

(259) 

(260) 

(261) 

(262) 


Perkins  and  Hage  (Reference  29)  present  a  formula  for  the  desirable  quantity  of 
for  satisfactory  handling  characteristics. 


(C„fl  )  3  0.0005 

'  "P  'desirod  ^  S 


57,  3  /  rod 


8 

(263) 
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When  the  last  two  equations  are  combined,  the  required  fin  area  may  be  ascertained 


I  I 


Sf  / 

’  'l*  1 

r  w  \ 

s  ■  °  ' 

^  s  ' 

(264) 


For  m  =  0.3  (Section  3-C)  and 


W 

S  IP- 


=  50  (Equation  253), 


0.1  . 


Thus,  the  fin  area  requirement  becomes  the  low  speed  stability  constraint. 


(265) 


77 


ASl)-'rDR-62-lI02 


LIST  OF  KCFiiRliNCES 


1.  Forsyth,  A.R.;  Calct'us  of  Variations:  Cambridge  University  Press,  1927. 

2.  Lees,  Lester;  "Laminar  Heat  Transfer  Over  Blunt  Nosed  Bodies  at  Hypersonic 
Flight  Speeds,'*  Jet  Propulsion;  Vol  26;  April  1956;  No.  4,  p.259. 

3.  Creager,  Marcus  O.;  Effects  of  Leading-Edge  Blunting  on  the  Local  Heat  Transfer 
and  Pressure  Distributions  Over  Flat  Plates  in  Supersonic  Flow;  NASA  TN  4142, 
National  Aeronautics  and  Space  Administration,  1957.  ~ 

4.  Wagner,  Pine,  and  Henderson;  Laminar  Heat  Transfer  and  Pressure  Distribution 
Studies  on  a  Series  of  Re-entry  Nose  Shapes  at  a  Mach  Number  of  19.4  in  Helium; 
NASA  TN  D-891i  National  Aeronautics  and  Space  Administration,  June  1961. 

5.  Hankoy,  W.L.,  Neumann,  R.D.,  Fllnn,  E.V.;  Design  Procedures  for  Computing  Aero¬ 
dynamic  Heating  at  Hypersonic  Speeds:  WADC  TR  59-610,  Wright  Air  Development 
Center,  Wright -Patterson  Air  Force  Base.  Ohio,  1960. 

6.  Eckert,  E.R.G.;  Survey  of  Heat  Transfer  at  High  Sneeds:  WADC  TR  54-70,  Wright 
Air  Development  Center,  Wright-Patterson  Air  Force  Base,  Ohio,  1954. 

7. ’  Schmidt,  J.F.:  Laminar  Skin  Friction  and  Heat  Transfer  Parameters  for  a  Flat 

Plate  at  Hypersonic  Speeds  in  Terms  of  Free-Stream  Flow  Properties:  NASA 
TN  D-8,  National  Aeronautics  and  Space  Administration,  September  1959. 

8.  Schmidt,  J.F.;  Turbulent  Skin  Friction  and  Heat  Transfer  Coefficients  for  an  Inclined 
Flat  Plate  at  High  Hypersonic  Speeds  in  Terms  of  Free-Stream  Flow  Properties; 
NASA  TN  D-869,  National  Aeronautics  and  Space  Administration,  May  1961. 

9.  Krier,  C.A.;  Coatings  for  the  Protection  of  Refractory  Metals  from  Oxidation:  DMIC 
Report  162;  Battalle  Memorial  Institute,  Columbus,  Ohio,  Noyember  1961. 

10.  Hayes,  Wallace  D,,  and  Probsteln,  R.F.;  Hypersonic  Flow  Theory;  Academii.  Press, 
New  York,  1959,  Chapter  111. 

11.  Sl  aniro,  Ascher  H.;  The  Dynamics  and  Thermodynamics  of  Compressible  Fluid  Flow; 
Volumes  1  and  11,  The  Ronald  Press  Company,  New  York,  1953. 

12.  Grey,  Don  J.:  Drag  and  Stability  Derivatives  of  Missile  Components  According  to 
Modified  Newtonian  Theory;  AEDC  TN-60  191,  November  1960. 

13.  Grimminger,  G.,  Williams,  E.P.,  and  Young,  G.B.S.;  "Lift  on  Inclined  Bodies  of 
Revolution  in  Hypersonic  Flow;”  Journal  of  Aeronautical  Sciences.  Volume  17,  No. 

11,  p.  675-690,  November  1950. 

14.  Lees,  Lester;  "Hypersonic  Flow,”  Proceedings  lAS-RAeS  Fifth  International  Aero¬ 
nautical  Conference,  Los  Angeles,  California.  June  20-23,  1955. 

15.  Lees,  L.;  "Recent  Developments  In  Hypersonic  Flow;”  Jet  Propulsion,  p.  1162- 
1178,  November  1957. 


79 


ASD-TDR-62-1102 


LIST  OF  RliFERFNJFS  (CCWT’D) 


16.  Love,  F.S.:  “Generalized  Newtonian  Theory;"  Journal  of  the  Aero/Space  Sciences, 
p.  314,  May  L)59. 

17.  Kaufman.  Louis  G.  11.  and  Scheuing,  Richard  A.;  “An  Introduction  to  Hypersonics," 
Grumman  Aircraft  Engineering  Corporation  Report,  RE-82,  October  1956. 

18.  Truitt,  Roliert  W.;  Hypersonic  Aerodynamics;  The  Ronald  Press  Company,  New 
York,  1959. 

19.  Malvcstuto,  Frank  S.,  Jr.,  et  al;  Study  to  Determine  Aerodynamic  Characteristics  on 
Hypersonic  Re-Entry  Configurations:  WADD  Technical  Report  61-56,  Part  II,  Volume 
r.  Wright  Air  Development  Division.  Wright-Patten-on  Air  Force  Base,  Ohio,  August 
1962. 

20.  Feldman,  Saul;  Hypersonic  Gas  Dynamic  Charts  for  Equilibrium  Air;  AVCO  Research 
Lab.,  AVCO  Manufacturing  Corporation,  January  1957. 

21.  Equations.  Tables,  and  Charts  for  Compressible  Flow;  NASA  Report  1135,  National 
Aeronautics  and  Space  Administration,  1953. 

22.  Manellucci,  A.,  and  Fields,  A.;  Pressure  Distribution  on  Blunted  Flat  Plates  with 
Surface  Inclination:  General  Applied  Science  Laboratories  Data  Report  No.  5134-5, 
Technical  Report  No.  171,  June  1960. 

23.  Creager,  M.O.:  Effects  of  Leading-Edge  Blunting  on  the  Local  Heat  Transfer  and 
Pressure  Distributions  Over  Flat  Plates  in  Supersonic  Flow;  NASA  TN  4142,  National 
Aeronautics  and  Space  Administration,  December  1957. 

24.  Creager,  M.O.;  The  Effect  of  Leading-Edge  Sweep  and  Surface  Inclination  on  the 
Hypersonic  Flow  Field  Over  a  Blunt  Flat  Plate;  NASA  Memo  12-25-58A,  National 
Aeronautics  and  Space  Administration,  January  1959. 

25.  Brown.  R.C..  et  al:  Six-Degree  of  Freedom  Flight  Path  Study  Generalized  Computer 
Program,  Part  1.  Problem  Formulation:  WADD  TP  60-781,  Wright  Air  Development 
Division,  Wright-Patterson  Air  Force  Base,  Ohio,  November  1960. 

26.  Etkin,  Bernard:  Dynamics  of  Flight,  John  Wiley  &  Sons,  Inc.,  New  York,  1959. 

27.  Matranga,  Gene  J.;  Analysis  of  X-15  Landing  Approach  and  Flare  Characteristics 
Determined  from  the  First  30  Flights.  NASA  TN  D-1057,  National  Aeronautics 
and  Space  Administration.  1961. 

28.  Donovan,  A.F.  (Editor);  Aerodynamic  Components  of  Aircraft  at  High  Speeds  Volume 
vn  fHlgh  Speed  Aerodynamics  and  Jet  Propulsion);  Princeton  University  Press, 

N.  J.,  1957. 

29.  Perkins,  C.D.,  and  Hage,  R.E.;  Airplane  Performance  StaPlllty  and  Control.  John 
Wiley  &  Sons,  Inc.,  New  York,  1949. 


80 


Tlie  complete  geometry  for  one  of  the  typical  optimums  The  complete  geometry  for  one  of  the  t>'pical  optimums 

U  shown  as  ijj  example*  shown  as  an  example. 


,1 

f 


iJ  c 

^  o  < 
fci 


<  <  5 


*3  H  o 

w  53  o  < 
g  =  i:ti 

H?--! 

i*  i  c  5 


a  3  -; 


imrsfimvtpPiMMni 


•V  H 

?  -  t 

d  ^  u 
u  ^  . 

3  c.  y 
a  §  S-  i 

V  u  S  c 

1  ^  5  t 

2  'i:  <  c2 


^  ^  -  a  3  > 


o 


o“  .Si  5  / 

1 1 

f£a  : 

•3  ^  tt.  ■* 
u  ^  a  3 

s  §  2  5 

.s 

“  a  <  . 

a  “r  i*  § 

^^i:- 

a  C  fc  3 

I  2 
|ir; 

®  •  «  VJ 
«  *5  I  u: 
c  J5  'sj  1 

ga-f 
■S  S  “  I 

llis 

£ 


o: 

t  • 
&«  « 


I  ?■?  «  a‘i  ?  .  -3  2 

|il|l|si'll 

n  =  i9li):y 


£  2  ' 
Ijy  V 


u  •■ajsig’S 

£|-3S2fSf9« 
“  3  I! 


ais 

e  - 


^  .s  ,s  £  Si  a 

5 I  Is 

^  ^  e  44  tt> 


-  1 1 
*  o  3 


A  ^  ■ 


HI  a  2  s I E -|3 

sl1l!|l5»| 

?•  «  S  5  J  g  S 


S  p  s 


-«Es  ii:s||1»li 

"'iwSEcigg.Ss^S.i 


i -SI  all  ill:  1 5 

s&Sf<io^55^cF2^<fl:S.SS 

5ui|ISsi?S|a|»Sp| 


t»  i' «  *0 

*1Si 

i  =  >  d 


■  in  c 

!  i  &l 
:  £  2  ? 

I  a  =  ='  -  •  > 

!  <  <  a  ^  s  <  3 


2  w  H 


'  =  3 : 


g 


^  O  O  9) 

a  2  5  j 

n.tt 

f£3  : 

a  ^  « 

"  I  °  i 

“£  fc*  O  o 

c  «{  C  £ 
“  a  5  . 
a'^'  2 1 
:itS* 
E 

8  it  o 


S| 

O 


*  O  i. 

>  >  2 

B  <M  % 


'iS 


•  *5  .  - . 

§  O  o 

tS  f  B  X 

y  5:^ 

vA^i 

imi 


sz 


£•  » 

.,  3  o 

4  t 


slllfligll,' 

r  »S  .=?«  I  S  S 

i?|lliii|f 

iltlilllll 

'■llllll 

lillllls 


s  t  s 

J5  ®  C  •  ,3 

ellill 

If S  8£| 


'  ^  h  •>  S  S 

|:l||yii 
t!  si 
ill  ill 


SaH  S*  S 

a » s  5 ,  -  -2  “•» 

I g i l|| g II 

.||  |i||£  . 

ii^la 


i«  y  ft! 

IlSIia; 


'V 


UNCLASSIFIED 


UNCLASSIFIED 


